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Abstract 


The  time  dependent  deformation  of  coal  under  a  constant 
load  can  be  described  adequately  by  the  sum  of  two  power 
laws : 

B  D 

e  =  At  +  Ct 

where  e  is  the  strain  rate 

t  is  the  elapsed  time  since  the  material  has 
been  under  the  constant  load 
A,  B ,  C  and  D  are  constants. 

The  attempt  to  find  the  stress  dependence  of  the  creep 
rate  of  a  model  coal  pillar  was  frustrated  by  the 
variability  of  the  material.  The  Young’s  Modulus,  E,  of  coal 
can  vary  widely  within  a  relatively  small  area,  which 
indicates  that  samples  from  the  same  sampling  site  can  have 
much  different  deformation  properties. 

Cruden  (1983)  suggested  that  the  stress  dependence  of 
creep  rate  can  be  determined  by  an  increment  test  on  a 
single  rock  specimen.  Such  a  test  eliminates  variation 
caused  by  variation  of  physical  properties  between 
specimens.  Further  research,  using  this  type  of  test  might 
allow  the  determination  of  the  physical  parameters 
controlling  the  creep  of  coal. 
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1 .  Introduction 


1.1  Introductory  Remarks 

The  study  of  time-dependent  effects,  usually  spoken  of 
under  the  general  title  of  'creep',  is  of  the  greatest 
importance  in  rock  mechanics  and  geophysics.  The  phenomenon 
of  time-dependent  behaviour  of  rocks  is  a  source  of  many 
problems  in  designing  structures  in  rock.  In  underground 
work,  the  movements  which  occur  after  excavation  are  of 
creep  type  and  the  requirement  is  to  find  laws  by  which 
future  behaviour  may  be  predicted. 

Since  the  original  work  by  Andrade  (1910),  who  studied 
the  behaviour  of  metal  wires  subjected  to  constant  tensile 
stress  above  the  elastic  limits,  there  has  been  developments 
of  the  time-dependent  theory  covering  both  a  wide  range  of 
materials  (such  as  metals,  rubber,  ice,  soils  and  rocks)  and 
applications  (such  as  civil,  mining  and  mechanical 
engineering,  etc.).  There  are  many  useful  reviews  of  creep 
of  rocks  by  Robertson  (1964)  and  Murrell  and  Misra  (1961) 
and  others.  Time-dependent  behaviour  in  several  rocks  has 
been  studied  both  in  the  laboratory  and  in-situ.  The 
decelerating  strain  rate  period  of  creep  has  been  well 
studied  and  several  creep  laws  were  developed,  but  there 
were  no  report  on  the  accelerating  strain  rate  period  known 
to  the  author. 

There  are  several  creep  equations  for  rocks  published 
in  for  example  Obert  and  Duvall  (1967),  Cruden  (1971), 
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Jaeger  (1972).  The  power  law  proposed  by  Cruden  (1970)  is 
used  in  this  thesis  for  two  reasons:  firstly,  there  are 
existing  programs  to  reduce  and  fit  the  experimental  data 
into  the  power  law;  and  secondly,  to  present  the  hypothesis 
that  accelerating  strain  rate  period  also  follow  the  power 
law.  Coal  was  chosen  as  the  testing  material  because  it  is  a 
rock-like  material. 

1.2  The  Objective  of  the  Study 

The  objectives  of  this  study  are: 

a.  To  determine  the  functional  form  of  the  creep  of  a 
model  coal  pillar. 

b.  To  confirm  the  functional  form  of  creep  by  a  large 
number  of  short  experiments. 

c.  To  use  these  experiments  to  check  whether  there  is  a 
critical  creep  strain  for  model  coal  pillars. 

1.3  Structure  of  Thesis 

Because  the  basic  creep  theory  can  be  found  in  many 
rock  mechanics  books,  therefore  it  is  not  presented  in  great 
detail  here.  Also,  no  detailed  review  on  developments  in 
creep  theories  is  included  as  it  had  been  done  in  previous 
research  work  by  da  Fontoura  (1980). 

Chapter  two  consists  of  a  description  of  the  testing 
program.  It  gives  a  detailed  account  of  the  sampling  site 
and  procedures.  The  method  used  in  sample  characterization 
and  the  testing  equipments  are  also  described. 
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Chapter  three  reports  the  method  used  in  data  analysis 
and  problems  that  had  been  encountered.  Results  from  the 
Star-Key  Coal  Mine  test,  which  was  conducted  by  Jeremic 
(Personal  Communication),  are  also  presented. 

In  Chapter  four,  results  from  the  testing  program  are 
presented.  A  comparison  of  test  results  to  those  reported  by 
da  Fontoura  is  also  included. 

Finally  in  Chapter  five,  discussion  of  test  results  and 
conclusions  are  presented. 


2.  Preparations  for  Laboratory  Tests 


2.1  Sampling  Site  and  Procedures 

The  coal  samples  used  in  the  present  study  were 
obtained  from  the  coal  seams  exploited  at  the  Highvale  Mine, 
which  is  situated  on  the  south  shore  of  Wabamun  Lake.  The 
Wabamun  Lake  district  is  west  of  Edmonton,  Alberta  in  Tps. 
50-54,  Rs.  3-7,  W.  5th  Mer.,  and  is  centered  about  Wabamun 
Lake . 

The  geology  as  well  as  the  topography  and  drainage  of 
the  area  in  the  proximity  of  the  sampling  site  have  been 
described  by  Pearson  (1959)  and  Noonan  (1972).  The  coal  unit 
is  referred  to  as  the  Pembina  Coal-bearing  zone  by  Pearson 
(1959),  and  is  sub-divided  into  several  seams  in  the 
vicinity  of  Wabamun  Lake.  The  coal-bearing  unit  can  be 
divided  into  two  main  seams  with  a  few  thinner  seams  below. 
The  two  main  seams  are  generally  about  three  metres  thick 
and  are  separated  by  an  interval,  from  a  few  centimetres  to 
ten  metres  of  shale  and  sandstone.  The  thick  seams  have  been 
termed  the  Upper  Main  and  Lower  Main  seams  (Pearson,  1959). 

The  coal  seams  at  Highvale  Mine  are  exploited  by  a 
conventional  strip-mining  operation.  The  till  cover  is 
removed  by  a  dragline  leaving  the  coal  seam  exposed  and 
light  explosive  charges  are  set  in  boreholes  to  loosen  the 
coal,  thereby  facilitating  the  mining  operation. 

Observations  of  the  blast  holes  exposed  along  the  face 
of  the  bench  were  made  by  Noonan  (1972),  who  indicated  that 
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the  visually  detectable  shatter-zone  extended  in  a  fan-like 
arrangement  only  about  45  cm  from  the  point  where  the  charge 
was  detonated.  Hence,  only  few,  if  any,  additional  fractures 
would  be  created  at  the  top  of  the  seam  as  a  result  of 
blasting . 

Sampling  was  carried  out  on  top  of  the  exposed  Lower 
Main  seam.  A  water-operated  laboratory  drilling  machine, 
manufactured  by  the  Milwaukee  Electric  Tool  Corp. ,  with  a 
core  barrel  of  about  7.5  cm  in  external  diameter  was  used  in 
the  coring  operations.  Reaction  against  the  weight  of  the 
field  vehicle  was  provided  for  the  drilling  machine  in  order 
to  avoid  unwanted  vibration  of  the  core  barrel  that  could 
damage  the  core.  Cylindrical  cores  of  about  7  cm  in  diameter 
were  drilled  in  the  site. 

Problems  were  encountered  during  the  coring  operations. 
The  fracturing  in  zones  inside  the  coal  seams  prevented  the 
successful  coring  of  samples  because  of  breakage  of  the 
core.  Also,  partings  inside  the  coal  seams  caused  cores  to 
separate  into  lengths  which  were  less  than  the  minimum 
requirement  of  length  to  diameter  ratio  of  2  to  1.  Only  one 
out  of  three  sampling  trips  was  successful  where  the  coal 
seam  is  more  intact  and  relatively  free  of  partings.  Two 
days  were  spent  for  the  successful  sampling  trip  and  a  total 
of  twenty-seven  samples  were  obtained. 
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2.2  Sample  Preparation 

After  the  samples  were  drilled  at  the  site,  the  cores 
were  removed  manually  by  pushing  them  from  the  core  barrel 
and  carefully  wrapped  with  Saran  Wrap  and  aluminum  foil  to 
prevent  moisture  loss.  Samples  were  then  put  into  boxes  and 
field  work  clothes  were  put  between  layers  of  samples  to 
minimize  the  possibility  of  breakage  during  transportation 
to  the  laboratory. 

After  transporting  the  test  samples  to  the  laboratory, 
they  were  stored  in  a  moist  room  at  5°C  and  100%  relative 
humidity  to  await  trimming.  The  samples  were  cut  into 
desired  lengths  using  a  Northland  concrete  saw  manufactured 
by  Oxford  Machine  and  Welding  Co.  Ltd.  of  Edmonton,  Alberta. 
The  saw  has  a  blade  diameter  of  60  cm  and  water  was  used  as 
cutting  fluid.  The  criterion  used  in  selecting  the  length  of 
the  sample  was  to  keep  the  length  to  diameter  ratio  at 
around  2  to  1 .  Special  wooden  holders  were  manufactured  to 
hold  the  samples  in  the  cutting  machine.  The  wooden  holders 
were  made  from  two  pieces  of  wood  of  size  200  mm  W  x  150  mm 
H  x  63.5  mm  T.  A  segment  of  circle,  with  the  same  diameter 
as  the  samples,  was  cut  from  the  face  of  each  piece  of  wood 
in  order  to  accomodate  the  sample.  Finally,  the  wooden 
holders  were  covered  with  a  layer  of  roofing  tar,  a 
water-proof  material,  to  prevent  the  swelling  of  the  wood. 
Coal  samples  were  first  wrapped  with  a  rubber  membrane 
before  being  put  in  between  the  wooden  holders.  The  rubber 
membrane,  the  same  type  used  for  any  triaxial  tests,  acted 
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as  a  shock-absorbing  medium  to  minimize  unwanted  vibrations 
while  samples  were  being  cut  with  the  saw  to  provide  some 
confinement  of  sample. 

After  fine-tuning  the  adjustable  alignment  of  the 
concrete  saw  and  the  sample  holders,  samples  could  be  cut 
with  acceptable  parallelism  of  the  end  surfaces.  The 
parallelism  was  measured  by  a  dial  gauge  which  was  connected 
to  a  smooth  level  platform.  By  turning  the  sample  on  the 
smooth  level  platform  with  the  dial  gauge  arm  touching  the 
sample’s  top  end,  the  maximum  difference  of  parallelism 
could  be  determined.  All  the  samples  used  in  testing  had 
their  ends  parallel  to  each  other  within  the  tolerance  limit 
of  0.25°,  i.e.,  to  within  approximately  5  mm/m  (Pit  Slope 
Manual,  Supplement  3-5,  1977). 

2.3  Sample  Characterization 

There  were  two  different  sources  of  coal  samples  used 
in  this  study.  The  first  source  was  samples  obtained  from 
the  sampling  site,  which  were  drilled  perpendicular  to  the 
bedding  planes,  as  illustrated  in  Figure  2.1.  The  second 
source  was  samples  not  used  by  da  Fontoura  (1980),  which 
were  drilled  from  a  block  sample  with  their  long  axis 
parallel  to  the  bedding  planes  and  at  an  angle  of  30°  to  the 
major  cleat,  as  shown  in  Figure  2.2.  A  detailed  structural 
survey,  conducted  by  Noonan  (1972)  concluded  that  the  major 
cleat  in  the  coal  seam  is  oriented  approximately  N45°E.  Thus 
the  configuration  of  the  samples  from  the  second  source 
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Figure  2.1  Coal  structure  of  samples  from  the  sampling  site 
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Figure  2.2  Coal  structure  of  da  Fontoura's  samples 
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would  correspond  to  a  sample  with  its  core  axis  horizontal 
and  oriented  either  at  N15°E  or  N75°E  in  the  field. 

Both  types  of  samples  were  prepared  and  cut  into 
desired  lengths  as  outlined  in  Section  2.2.  Specimen  weight 
and  dimensions  were  recorded,  and  sample  volume  and  density 
were  then  calculated.  The  samples  were  characterized  by 
measuring  the  velocities  of  compressional ,  Vp,  and  shear. 
Vs,  elastic  waves  through  each  sample.  The  Young's  Modulus, 
E,  of  each  sample  could  then  be  determined  from  the 
following  equations  (Jaeger,  1962): 


Vp=[ (X+2G)/p] 1 ' 2 
Vs= (G/p ) 1 7  2 
E=G ( 3X+2G ) /( X+G ) 


(2.1  ) 


(2.2) 


(2.3) 


where  X  and  G  are  Lame's  parameters  and  p  is  the  density  of 
the  sample. 

The  velocities  of  compressional,  Vp,  and  shear,  Vs, 
elastic  waves  were  determined  by  dividing  the  length  of  the 
sample  by  the  time  for  the  respective  wave  to  travel  through 
the  sample.  Each  sample  was  clamped  between  a  set  of 
compressional  sonic  heads  and  shear  sonic  heads  in  order  to 
measure  the  time  for  the  compressional  and  shear  elastic 
waves,  respectively,  to  travel  through  the  sample.  Vaseline 
was  used  as  a  conducting  agent  between  the  heads  and  the 
sample  in  both  cases.  Sonic  waves  were  sent  from  one  head  to 
the  other  through  the  length  of  the  sample  by  a  Terrametrics 
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Sonic  Pulse  Generator.  The  times  were  measured  by  a 
oscilloscope  manufactured  by  Gould  Advance  Ltd. ,  England. 

Before  the  apparatus  were  used  to  determine  the 
velocities  of  compressional  and  shear  elastic  waves  through 
the  coal  samples,  a  calibration  test  was  done  by  using 
aluminum  cores.  Aluminum  cores  of  five,  ten  and  twenty 
centimetres  in  length  were  used  and  the  time  of 
compressional  and  shear  wave  travel  were  determined  in  the 
same  way  as  outlined  above.  Vaseline  was  also  used  as 
conducting  agent  between  the  heads  and  the  core.  Plots  of 
distance  travelled  against  time  were  made  for  compressional 
and  shear  waves,  and  the  velocities  were  determined  as  the 
slope  of  the  straight  line  of  the  respective  plots.  The 
Young's  Modulus,  E,  of  aluminum  was  determined  by  using 
Equations  (2.1)  to  (2.3).  The  calculated  Young's  Modulus  of 
aluminum  is  within  4%  of  the  published  value. 

Despite  the  accuracy  in  determining  the  Young's  Modulus 
of  aluminum,  errors  could  arise  in  many  different  ways.  The 
time  of  travel  of  waves  is  slightly  affected  by  the  amount 
of  vaseline  used  and  also  the  clamping  force.  The  time  of 
travel  is  also  subject  to  error  in  reading  the  time  from  the 
oscilloscope.  Since  coal  is  a  heterogeneous  material  due  to 
its  composition  and  because  of  its  discontinuous  nature 
(Kaiser  and  Maloney,  1982),  it  is  more  difficult  to 
determine  the  time  of  travel  through  coal  than  aluminum. 

Time  error  accounts  for  most  of  the  errors  that  arise  in 
sample  characterization.  From  Equations  (2.1)  to  (2.3),  the 
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Young's  Modulus,  E,  can  be  written  in  terms  of  Vp  and  Vs  as 


E  =  pV s  2 ( 3Vp2  ~4Vs 2 )  (2.4) 

Vp2-Vs2 

It  is  obvious  that  the  value  of  Vs  can  affect  the  value  of 
the  Young's  Modulus,  E,  to  a  great  extent.  Also,  the  error 
in  Young's  Modulus  is  greater  than  the  error  in  Vs  as  E  is 
directly  proportional  to  the  square  of  Vs. 

The  results  of  sample  characterization  are  summarized 
in  Table  2.1.  The  last  five  samples  in  Table  2.1,  T6A  to 
T6J,  are  samples  left  over  from  da  Fontoura  (1980).  Since 
these  five  samples  were  drilled  in  a  different  orientation 
from  the  others,  and  the  differences  were  reflected  in  the 
results  of  sample  characterization.  Da  Fontoura' s  samples 
were  drilled  parallel  to  the  bedding  planes  and  yielded 
higher  values  of  compressional  and  shear  wave  velocities  as 
the  waves  propagate  parallel  to  the  bedding  and  only  have  to 
cross  the  joints.  Samples  from  the  sampling  site  were 
drilled  perpendicular  to  the  bedding  planes,  they  yielded 
lower  values  of  compressional  and  shear  wave  velocities 
because  the  waves  have  to  propagate  across  the  bedding 
planes.  Samples  drilled  parallel  to  the  bedding  planes  tend 
to  fail  in  shear  failure  along  the  joints  or  by  buckling 
along  individual  bedding  planes.  For  samples  drilled 
perpendicular  to  bedding  planes,  cracks  will  be  closed  by 
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SAMPLE  HEIGHT  DENSITY  Vp  Vs  E 

NUMBER  (mm)  (gm/cm3)  (mm/ sec)*10**6  (MPa) 
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Table  2.1  -  Summary  of  Sample  Characterization 


r* 


J  -•  O  O  -  —  — 


'  •  • 


14 


compression  and  rupture  by  crushing. 

2.4  Testing  Equipment 

A  simple  double-lever  arm  rig  capable  of  applying  a 
constant  axial  load,  maximum  capacity  of  183  kN,  was  used 
for  the  series  of  creep  tests  reported  here.  The 
double-lever  arm  rig  outlined  in  Figure  2.3,  was  originally 
designed  and  built  for  da  Fontoura  (1980).  The  rig  consists 
of  a  reaction  frame  and  two  lever-arms  (I-Section)  which 
would  transfer  loads  applied  at  their  ends  through  a  loading 
ram  to  the  sample.  The  mechanical  magnification  for  the 
double-lever  arm  system  was  7.5.  Before  loading  the  sample, 
hydraulic  jacks  were  used  to  support  the  weight  at  the  ends 
of  the  lever-arms.  The  positions  of  the  two  hydraulic  jacks 
is  illustrated  in  Figure  2.1.  Sudden  loading  was  achieved  by 
releasing  the  hydraulic  jacks  simultaneously. 

The  creep  rig  presented  a  problem  when  high  axial  load 
was  required,  large  number  of  weights  had  to  be  put  on  the 
hangers.  In  order  to  achieve  sudden  loading,  the  hydraulic 
jacks  supporting  the  weights  had  to  be  released 
simultaneously.  If  the  hydraulic  jacks  were  released  one  at 
a  time,  the  weights  on  one  of  the  hangers  would  cause 
instability  of  the  creep  rig.  It  would  be  desirable  to  widen 
or  stabilize  the  base  of  the  creep  rig  and  thus  to  eliminate 
the  instability  problem. 

A  triaxial  cell  for  10  cm  diameter  samples  was  modified 
in  order  to  accomodate  7  cm  diameter  samples  by  changing 
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Figure  2.3  Sketch  of  creep  rig 
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both  the  top  cap  and  the  bottom  pedestal.  Special  Thompson 
linear  bushings  were  used  to  guide  the  loading  ram  with 
minimal  shaft  friction.  The  triaxial  cell  used  had  drainage 
for  the  sample  provided  at  both  top  cap  and  bottom  pedestal. 

A  unit  for  monitoring  axial  load,  displacement,  and 
confining  pressure  complemented  the  laboratory  set-up.  This 
unit  is  a  Hewlett  Packard  3054  data  logger  which  consists  of 
a  HP3497A  data  aquisition  unit  and  a  HP85  computer.  The  HP85 
computer,  which  can  be  programed  to  take  readings  at  a 
preset  time-interval,  has  a  built-in  tape  recording  device 
and  is  capable  of  storing  all  the  information  in  the  tape.  A 
more  efficient  HP82901M  Flexible  Disc  Drive  Unit  was  later 
introduced  to  replace  the  tape  recording  device.  The  HP85 
computer  is  compatible  with  the  Michigan  Terminal  System 
(MTS)  at  the  University  of  Alberta,  and  data  can  be 
transferred  directly  from  the  tape  or  disc  through  the  HP85 
computer  to  the  MTS  for  analysis. 

The  axial  displacements  of  the  samples  were  measured 
with  two  Linearly  Variable  Differential  Transformers  ( LVDTs ) 
manufactured  by  Hewlett  Packard.  Initially  two  24DCDT-050 
LVDTs  with  a  displacement  range  of  ±1.27  mm  were  used. 
However,  it  was  found  that  the  displacement  of  the  sample 
was  greater  than  the  displacement  range  of  the  LVDTs.  The 
two  DCDT-050  LVDTs  were  then  replaced  by  a  pair  of 
24DCDT-100  LVDTS  with  a  displacement  range  of  ±2.54  mm.  The 
axial  load  was  measured  with  a  Transducer  Load  Cell, 
manufactured  by  Transducer  Inc.,  California,  U.S.A.,  with  a 
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capacity  of  45  kN  (10,000  pounds).  The  confining  pressure 
was  measured  with  a  transducer,  manufactured  by  Celesco 
Trans  Product  Inc.,  California,  USA,  Model  PLC  with  a 
capacity  of  691  kPa  (100  psi).  A  power  supply  unit  capable 
of  providing  input  voltage  of  6  Volts  was  used  to  feed  the 
load  cell  and  the  transducer.  Another  power  supply  unit  was 
used  to  supply  a  24  volts  input  voltage  to  the  LVDTs . 

The  LVDTs,  load  cell  and  transducer  were  calibrated 
before  being  used  and  no  change  in  the  calibration  factors 
was  observed  during  the  experimental  program.  The  testing 
apparatus  was  kept  in  a  temperature  and  humidity-controlled 
room,  Room  B21A  in  the  Civil  Engineering  Building,  capable 
in  keeping  the  temperature  variations  within  2°C  and  the 
humidity  within  5%. 

Even  though  the  testing  room  was  temperature  and 
humidity  controlled,  there  were  slight  variation  in 
temperature  and  humidity  between  day  and  night.  Equipment 
with  high  sensitivity  such  as  the  LVDTs,  were  affected  the 
most  and  this  caused  the  scattering  of  data.  Generally,  from 
the  read  out  of  the  LVDTs,  the  deformation  of  a  sample  under 
stress  would  be  larger  during  the  day  and  smaller  during  the 
night.  This  problem  could  be  solved  with  a  more  effective 
and  sophisticated  way  of  temperature  and  humidity  control, 
however,  at  a  larger  expense. 
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2.5  Testing  Procedures 

There  were  three  different  types  of  tests  carried  out 
in  this  experimental  program.  The  first  type  was  tests  with 
confining  and  back  pressure,  referred  to  as  the  C-tests. 
Prior  to  set-up,  each  sample  was  enclosed  within  a  double 
rubber  membrane  as  an  extra  precaution  to  avoid  leakage  in 
case  one  membrane  was  punctured  during  the  test.  Double 
O-rings  and  screw  clamps  were  used  to  provide  extra  seals 
along  the  contacts  between  membrane  and  both  top  cap  and 
pedestal.  Both  confining  and  back  pressure  were  applied 
before  the  sample  was  loaded. 

The  second  type  of  test  was  the  uniaxial  compression 
test.  The  tests  were  carried  out  in  the  modified  triaxial 
cell  with  no  rubber  membrane  and  no  cell  fluid.  This  type  of 
tests  was  referred  to  as  the  U-tests. 

The  third  type  of  test  was  with  no  rubber  membrane  on 
the  sample  and  using  small  confining  pressure,  referred  to 
as  the  S-tests.  Prior  to  set-up,  each  sample  was  saturated 
by  submerging  it  in  water  for  about  24  hours.  This  type  of 
test  was  carried  out  because  the  samples  were  too  strong  for 
failure  under  short  term  creep,  thus  samples  had  to  be 
weakened  by  saturation. 

For  all  the  tests  reported  herein,  sudden  loading  was 
obtained  by  opening  the  valves  on  the  two  hydraulic  jacks 
that  were  supporting  the  weights.  The  duration  of  the  creep 
tests  ranged  from  three  days  to  more  than  two  weeks.  None  of 
the  samples  failed  under  the  constant  axial  load  during 
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creep  tests. 

In  the  testing  program,  some  samples  were  used  for 
several  creep  tests.  After  each  test,  the  sample  was 
carefully  removed  from  the  testing  apparatus  and  wrapped  in 
Saran  Wrap  and  aluminum  foil  to  prevent  moisture  loss.  A 
minimum  recovery  period  of  24  hours  was  allowed  for  samples 
to  be  used  in  successive  tests.  Thus,  each  test  performed  in 
the  testing  program  could  be  treated  as  individual  single 
stage  creep  test.  A  single  sample  was  used  for  tests  CT6A1 
to  CT6A9.  After  each  test  was  terminated,  the  sample  was 
unloaded  but  remained  in  the  triaxial  cell  for  24  hours 
before  another  test  was  carried  out  with  a  higher  load. 

Table  2.2  is  a  table  of  tests  that  were  carried  out  in 
this  testing  program,  the  load  and  duration  of  each  test  are 
also  indicated.  The  first  letter  of  the  name  of  each  test 
represents  the  type  of  test.  It  is  followed  by  the  sample 
number,  as  listed  in  Table  2.1.  The  number  following  the 
sample  number  represents  the  number  of  times  the  sample  is 
loaded.  For  example:  CT6A  is  a  confined  test  with  sample  T6A 
and  it  is  the  first  time  loading.  UT6F2  is  a  uniaxial 
compression  test  with  sample  T6F  and  it  is  the  second  time 
loading . 
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Table  2 


2  Table  of  tests 


Test 

Axial  Load  (kN) 

Dura  t 1  on 

CT  6  A  1 

7 

90 

CT6A2 

10 

12  1 

CT6A3 

10 

144 

CT6A4 

12 

120 

CT6A5 

15 

150 

CT6A6 

15 

14  1 

CT6A7 

15 

98 

CT6A8 

18 

67 

CT6A9 

2  1 

96 

C9-9*8 

15 

93 

C9-9*9 

20 

7  1 

U9-9*9 

25 

168 

U9-9*  1  1 

24 

70 

U9-9*  1  7 

16 

7  1 

U9-9*2  1 

24 

94 

U 10-9*2 

27 

92 

UT6F 

7 

70 

UT6C 

1  1 

70 

UT6D 

10 

72 

UT6F2 

9 

187 

UT6D2 

10 

7  1 

UT6D3 

1  1 

163 

59-9*7 

20 

360 

59-9*9 

17 

144 

S9-9*  1  1 

22 

142 

S9-9*  1  5 

15 

166 

S9-9*  17 

15 

166 

S9-9*2 1 

20 

142 

S 10-9*2 

16 

148 

( Hour  s ) 
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3.  Creep  Behaviour  from  Laboratory  Tests 


3.1  Analysis  of  Creep  Data 

There  are  two  steps  in  the  analysis  of  laboratory  tests 

which  are  of  equal  importance.  First,  the  presentation  and 

conditioning  of  the  experimental  observations  and  second, 

the  analysis  of  the  processed  data. 

The  displacement  of  the  sample  was  measured  at  a  number 

of  times  after  the  application  of  the  load,  it  was  then 

transformed  into  engineering  axial  strain,  e  ,  by  the 

1 

following  expression: 


e  =  (L  -  L  )/L  (3.1) 

•  • 

l  l 

where  L  represents  the  initial  length  of  the  sample  and  L 

i 

is  the  length  of  the  sample  at  the  time  when  the  reading  is 
taken . 

The  interpretation  of  creep  data  is  done  basically  in 
terms  of  strain  rate,  e,  which  is  the  change  of  total 
strain,  e,  per  unit  of  time.  The  definition  of  strain  rate 
can  be  expressed  mathematically  as  in  equation  (3.2), 

e  =  Ae/At  (3.2) 


Since  the  total  strain  is  known  only  at  certain  times, 

t  ,  the  estimation  of  the  strain  rate  has  to  be  done  by 
i 

numerical  differentiation.  The  simplest  approach  would  be  to 
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approximate  the  strain  rate,  e  ,  at  time  T  =  (t  +  t  )/2 

1  1  1  ~  1  1 

by  ( e  -  e  )/(t  -  t  ).  This  approach,  however,  presents 

i  i- 1  i  i- 1 

some  difficulties.  Small  fluctuations  in  the  output  voltage 

of  the  LVDTs  and  also  temperature  caused  some  observations 

of  strain,  e  ,  at  time  t  to  be  smaller  than  the 

•  • 

l  l 

observations,  e  ,  at  time  t  ,  which  corresponds  to  a 

i-1  i-1 

negative  strain  rate. 

Cruden  (1969)  proposed  to  smooth  the  original 

observations  using  recursion  formulae.  If  e  is  less  than 

i 

e  ,  a  new  observation  E  =  (e  +  e  ) /2  is  defined 

%  A  w  A  *  *  A 

1-1  1-1  1  1-1 

associated  with  a  time  T  =  (t  +  t  )/2.  The  new 

i-1  i  i-1 

observation,  E  ,  is  given  a  weight,  W  ,  which  is  equal 

i-1  i-1 

to  the  sum  of  w  and  w  .  For  the  original  data,  all 

i  i-1 

observations  have  a  weight,  w  ,  equal  to  unity.  This  process 

i 

is  followed  until  all  the  observations,  E  ,  are  such  that 

i 

every  strain  is  greater  than  the  previous  ones.  From  the  new 

set  of  observations,  (E  ,  T  ) ,  the  strain  rates  are 

i  i 

calculated  using  the  simple  approach  mentioned  earlier  in 
this  section. 

The  process  of  creep  deformation  can  be  divided  into 
two  main  regions.  Initially,  it  is  characterized  by  a 
decreasing  rate  of  strain,  this  is  called  the  decelerating 
creep  region.  There  follows  a  stage  where  the  rate  of  creep 
strains  increases  with  time  and  eventually  leading  to 
failure.  This  region  is  known  as  the  accelerating  creep 


region . 
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Cruden  (1971a)  suggested  that  the  decelerating  creep 
could  be  adequately  described  by  a  power  law.  The  power  law 
relationship  between  strain  rate  and  time  was  represented  by 
Equation  (3.3),  where  e  is  the  strain  rate,  A  is  a  constant 
which  is  stress  dependent  and  B  is  a  strain-hardening 
parameter . 

-B 

e  *  At  (3.3) 

This  equation  is  represented  by  a  straight  line  with  a 
negative  slope  in  a  double  logarithm  plot  of  strain  rate 
versus  time. 

It  is  believed  that  the  accelerating  creep  can  be 
represented  by  a  similar  power  law  except  with  a  positive 
power,  i.e.,  the  straight  line  representing  the  accelerating 
creep  in  a  double  logarithm  plot  of  strain  rate  versus  time 
will  have  a  positive  slope.  Thus,  the  entire  creep  curve 
might  be  represented  by  two  power  laws  as  in  equation  (3.4): 

3  D 

e  =  At  +  Ct  (3.4) 

where  the  first  term  on  the  right  hand  side  would  describe 
the  decelerating  creep  and  the  second  term  describes  the 
accelerating  creep. 
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3.2  Problems  Encountered  in  Data  Analysis 

Experimental  observations  obtained  from  laboratory 
testing  were  collected  and  processed  as  outlined  in  Section 
3.1.  The  time  and  strain  rate  were  calculated  accordingly. 

To  fit  the  experimental  data  into  the  form  of  Equation 
(3.4),  the  BMDP  Statistical  Software  was  employed.  The  BMDP 
Statistical  Software  is  prepared  by  the  Department  of 
Biomathematics,  University  of  California,  Los  Angeles.  It  is 
available  as  one  of  the  statistical  package  at  the 
University  of  Alberta. 

The  BMDP  computer  programs  (Dixon,  1981)  are  designed 
to  aid  data  analysis  by  providing  methods  ranging  from 
simple  data  display  and  description  to  advance  statistical 
techniques.  Data  are  usually  analyzed  by  an  iterative 
'examine  and  modify'  series  of  steps. 

There  are  two  BMDP  computer  programs  for  non-linear 
regression,  namely  the  P3R  and  PAR  programs.  To  use  these 
programs,  a  main  program  using  the  BMDP  Instruction  Language 
had  to  be  written.  The  main  program  gives  the  location  of 
the  input  data,  the  initial  estimates  of  the  parameters  A, 

B,  C  and  D  and  most  importantly,  the  specific  function 
required.  Table  3.1  gives  a  listing  of  a  sample  program 
using  the  BMDP  Instruction  Language. 

Program  P3R  gives  least  squares  estimates  of  the 
parameters  of  a  non-linear  function.  Six  functions  (and 
their  derivatives)  are  built-in.  Other  functions  can  be 
fitted  to  the  data  by  specifying  both  the  function  and  its 
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Table  3.1  Listing  of  a  sample  program  using  the  BMDP 
instruction  language 


/PROBLEM 

TITLE  IS  'STAR-KEY  COAL  MINE  NO. =5' 

/INPUT 

ERRLEV  IS  STRICT. 

VARIABLES  ARE  2. 

FORMAT  IS  FREE. 

/VARIABLE 

/REGRESS 

MTSFILE  IS  BMDPINPUT. 

NAMES  ARE  LNT ,  STRATE. 

DEPENDENT  IS  STRATE. 

/PARAMETER 

PARAMETERS  ARE  4. 

INITIAL  ARE  10.0,  -0.1,  0.5,  1.0. 

/  SAVE 

MINIMUM  ARE  0.0,  -1.0,  0.0001,  0.5. 
MTSFILE  IS  BMDPOUT PUT . 

CODE  IS  ' STARKEY5' . 

NEW. 

/END 
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derivatives  by  BMDP  Control  Language  statements  or  by 
FORTRAN  statements.  The  parameters  are  estimated  by  a 
Gauss-Newton  algorithm.  Upper  and  lower  limits  can  be  placed 
on  the  parameters,  and  exact  linear  constraints  for  the 
parameters  are  available. 

Similarly  to  program  P3R,  program  PAR  also  estimates 
the  parameters  of  a  non-linear  function  by  least  squares. 

The  program  is  appropriate  for  a  wide  variety  of  functions 
for  which  derivatives  are  difficult  to  specify  or  costly  to 
compute.  The  regression  function  must  be  specified  by 
FORTRAN  statements;  the  derivatives  are  not  specified.  Upper 
and  lower  limits  may  be  specified  on  the  individual 
parameters  or  for  arbitrary  linear  combinations  of  the 
parameters . 

In  order  to  use  either  program  P3R  or  program  PAR, 
Equation  (3.4)  had  to  be  transformed.  By  utilizing  the 
mathematical  equivalence  of  Equation  (3.5), 

B 

t  =  exp  (B  In  t)  (3.5) 

Equation  (3.4)  can  be  re-written  as 

e  =  A  exp  (B  In  t)  +  C  exp  (D  In  t)  (3.6) 

where  e  is  the  strain  rate 

A,  B,  C  and  D  are  constants 


♦ 
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Equation  (3.6)  becomes  the  fitted  form  of  the  power  law  and 
is  one  of  the  six  built-in  functions  in  program  P3R . 

Data  of  a  creep  test  for  the  Star-Key  Coal  Mine 
(Jeremic,  personal  communication)  were  analyzed  using 
programs  P3R  and  PAR  of  the  BMDP  Statistical  Software.  The 
experimental  observations  were  reduced  with  elapsed  times, 
the  time  difference  between  the  start  of  the  experiment  and 
the  time  that  the  observation  was  taken,  and  strain  rates 
calculated  as  outlined  in  Section  3.1.  Estimates  of 
parameters  A,  B,  C  and  D  were  obtained  by  separating  the 
decelerating  and  accelerating  creep  portions,  which  was 
achieved  by  fitting  a  best  fit  straight  line  to  the 
beginning  and  ending  portions  of  the  data  individually. 

These  estimates  were  used  as  initial  estimates  of  parameters 
in  progress  P3R  and  PAR.  Input  data  were  strain  rates  and 
the  natural  logarithm  of  time. 

Initially,  the  programs  were  run  with  no  upper  or  lower 
restraints  on  the  parameters.  A  convergence  problem  caused 
an  error  of  overflow  of  exponents  in  the  least  squares 
calculation  routine.  Different  measures  were  taken  to  solve 
the  problem,  including  using  double  precision  in  the  BMDP 
programs  and  setting  lower  limits  to  the  four  parameters. 

The  programs  were  finally  run  successfully,  but  the  results 
were  less  than  satisfactory.  Parameter  C  tends  to  approach 
zero,  which  makes  the  second  term  of  Equation  (3.6)  vanish. 

Programs  P3R  and  PAR  are  typical  of  existing  software 
programs  available  for  non-linear  regression.  It  is 
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concluded  that  the  programs  P3R  and  PAR  are  not  suitable  for 
analyzing  creep  data  because  these  programs  only  work  with 
data  with  small  scatters  (M.L.  Marshal,  Computing  Services, 
personal  communication).  Even  with  smoothing  techniques  for 
the  creep  data,  the  scatter  is  apparently  still  too  much  for 
the  BMDP  programs.  Thus,  another  method  of  data  analysis  had 
to  be  found. 

3.3  Method  of  Data  Analysis 

Since  the  programs  P3R  and  PAR  of  the  BMDP  Statistical 
Software  are  not  suitable  for  fitting  the  experimental  creep 
data  to  the  power  law  as  in  Equation  (3.4),  another  method 
which  involves  least  squares  regression  and  integration  of 
strain  rates  is  used. 

The  experimental  observations  were  used  to  calculate 
the  elapsed  time  and  strain  rates.  After  smoothing  the  data, 
they  were  fitted  into  a  decelerating  creep  power  law  by  the 
least  squares  method.  Then  the  decelerating  creep  power  law 
was  integrated  to  obtain  estimates  of  the  decelerating  creep 
strains.  The  accelerating  creep  strains  were  calculated  by 
subtracting  the  decelerating  creep  strains  from  the  observed 
strains.  From  the  accelerating  creep  strains,  accelerating 
strain  rates  were  calculated  and  fitted  to  an  accelerating 
creep  power  law  by  the  same  least  squares  method  used 
throughout  the  analysis.  As  we  will  see  from  the  parameters 
fitted  to  typical  experiments,  at  low  value  of  t  the 
contribution  of  the  accelerating  creep  curve  to  creep  is 
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negligible.  At  high  value  of  t,  the  contribution  of  the 
decelerating  creep  curve  to  creep  is  also  negligible. 

A  series  of  computer  programs  were  developed  to  carry 
out  the  analysis  as  outlined  above.  There  are  a  total  of  six 
computer  programs,  which  are  put  together  in  a  package, 
named  CPACK,  to  facilitate  data  analysis.  Programs  in  CPACK 
are  interactive  programs,  which  require  user  response  to 
terminal  prompts. 

The  first  program  in  package  CPACK  is  the  main  program 
which  outlines  the  steps  of  the  analysis  by  calling  the  five 
subroutines  in  the  package.  The  flow  diagram  of  CPACK  is 
shown  in  Figure  3.1.  The  five  subroutines  are  CRED,  BFIT, 
INTEG,  CFIT  and  CPLOT. 

The  main  program  first  called  the  Subroutine  CRED  to 
reduce  the  experimental  data,  compute  engineering  strains 
and  smooth  the  data  to  avoid  negative  strains.  The  flow 
diagram  of  Subroutine  CRED  is  shown  in  Figure  3.2. 

The  Subroutine  BFIT  is  then  called  to  convert  the 
strains  into  strain  rates,  and  to  fit  a  best  fit  straight 
line  to  either  the  decelerating  or  accelerating  creep 
portion  based  on  the  least  squares  method.  The  least  squares 
criteria  used  in  this  Subroutine  are  the  Durbin-Watson 
statistic  (Durbin  and  Watson,  1951)  and  the  Test  of  Slope 
Significance  (Cruden,  1971a).  These  two  criteria  are  to  be 
explained  in  more  detail  in  Section  3.4.  In  the  decelerating 
creep  portion,  the  reduced  data  are  fitted  by  a  straight 
line  with  negative  slope  in  a  double  logarithm  plot  of 
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Figure  3.1  Flow  diagram  of  program  CPACK 
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Figure  3 


2  Flow  diagram  of  Subroutine  CRED 
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strain  rate  versus  time.  When  the  least  squares  criteria  are 
not  satisfied,  the  last  data  point  is  taken  out  and  the  rest 
of  the  data  are  fitted  again.  The  process  is  repeated  until 
the  least  squares  criteria  are  satisfied.  The  flow  diagram 
of  Subroutine  BFIT  is  presented  in  Figure  3.3.  The  first 
call  to  Subroutine  BFIT  is  to  fit  a  best  fit  straight  line 
to  the  decelerating  creep  portion  in  a  double  logarithm  plot 
of  strain  rate  versus  time. 

The  best  fit  decelerating  strains  are  computed  in  the 
Subroutine  INTEG  by  integrating  the  power  law  expression 
fitted  for  the  decelerating  strain  rates  in  Subroutine  BFIT. 
The  decelerating  creep  strains  are  then  subtracted  from  the 
observed  creep  strains  to  obtain  the  accelerating  creep 
strains.  The  flow  diagram  of  Subroutine  INTEG  is  illustrated 
in  Figure  3.4. 

Subroutine  BFIT  is  called  again  to  convert  the 
accelerating  strains  computed  in  Subroutine  INTEG  into 
strain  rates  and  to  fit  a  power  law  to  the  accelerating 
creep  portion.  For  the  accelerating  creep  portion,  when  the 
least  squares  criteria  are  not  satified,  the  first  data 
point  is  taken  out  and  the  rest  of  the  data  are  fitted 
again.  This  process  is  repeated  until  the  least  squares 
criteria  are  satisfied. 

The  package  allows  the  user  to  have  the  option  of 
calling  Subroutine  CFIT,  which  calculates  the  strain  rates 
for  the  overall  fit,  by  combining  the  fits  of  the 
decelerating  and  accelerating  components.  It  then  calculates 
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Figure  3.3  Flow  diagram  of  Subroutine  BFIT 
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Figure  3.4  Flow  diagram  Of  Subroutine  INTEG 
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the  ratio,  R,  of  the  squares  of  the  scatter  of  the  data 
points  about  their  mean  to  the  squares  of  the  scatter  about 
this  overall  fit,  a  measure  of  goodness  of  fit.  Control  is 
then  transferred  back  to  Subroutine  EFIT  to  find  the  next 
smallest  range  of  data  which  satisfies  the  least  square 
criteria  with  a  new  power  law.  The  resulting  accelerating 
fit  parameters  are  transferred  to  Subroutine  CFIT,  which 
evaluates  the  overall  fit  obtained  with  the  new  parameters. 
The  process  is  repeated  until  the  overlap  between  the 
accelerating  and  decelerating  creep  goes  to  zero,  or  until 
the  accelerating  creep  strain  rate  comes  to  within  one  third 
of  the  decelerating  strain  rate  at  the  beginning  of  the 
range  of  data  used  for  the  fit  of  accelerating  creep.  The 
parameters  and  statistics  corresponding  to  each  trial  fit 
are  tabulated,  to  allow  the  user  to  identify  the  best  fit. 
The  flow  diagram  of  Subroutine  CFIT  is  shown  in  Figure  3.5. 

Finally,  Subroutine  CPLOT  is  called  to  produce  two 
plots  of  the  data:  time  versus  strain  rate  and  log  (time) 
versus  log  (strain  rate).  The  user  has  the  option  of  calling 
Subroutine  CPLOT  to  plot  the  data  and  the  best  fit  straight 
line  for  only  the  decelerating  creep  portion  or  the  overall 
fit  in  the  double  logarithm  plot  of  strain  rate  versus  time. 
The  flow  diagrams  of  Subroutine  CPLOT  are  illustrated  in 
Figure  3.6  and  3.7. 

This  method  of  data  analysis  works  very  well  for  data 
containing  both  decelerating  and  accelerating  creep  rate 
portions.  For  data  that  do  not  have  any  accelerating  creep 
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Figure  3.5  Flow  diagram  of  Subroutine  CFIT 
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Figure  3.6  Flow  diagram  of  Subroutine  CPLOT 
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Figure  3.7  Flow  diagram  of  Subroutine  CPLOT  (con't) 
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rate  portion,  an  error  message  will  likely  be  encountered  at 
the  second  call  of  Subroutine  BFIT .  The  cause  of  the  error 
may  be  due  to  too  few  points  left  after  the  decelerating  fit 
or  the  points  left  are  too  scattered  to  get  a  good 
accelerating  fit  that  satisfies  the  least  squares  criteria. 

The  user's  manual  of  the  computer  package  CPACK 
together  with  the  program  listings  are  put  together  in 
Appendix  A.  In  Section  3.4,  the  program  is  used  to  analyze  a 
long  (159  days)  creep  experiment. 

3.4  The  Creep  of  Star-Key  Coal 

The  Star-Key  Coal  Mine  is  located  at  Lsd.  4,  Sec.  36, 
Tp.  54,  R.  25,  W.  4th  Mer.,  approximately  17.7  km  north  of 
Edmonton,  Alberta  on  the  west  bank  of  the  Sturgeon  River 
valley.  The  coal  is  sub-bituminous  (Jeremic  and  Cruden, 
1979),  the  same  classification  as  the  coal  from  the  Wabamun 
Lake  District  (Pearson,  1959). 

Samples  were  brought  in  drums  from  the  mine  site  to  a 
laboratory  in  the  Mineral  Engineering  Department  at  the 
University  of  Alberta  to  protect  them  from  disintegrating 
during  transportation.  The  test  data  were  obtained  from  a 
particular  creep  test  carried  out  by  Jeremic  (Personal 
Communication).  The  test  was  a  uniaxial  compression  creep 
test  on  a  right  rectangular  prism  of  46.2  Wide  x  47.8  Long  x 
48.5  Height,  all  measurements  are  in  millimetres.  The  axial 
stress  was  reported  as  19.4  MPa  and  the  duration  of  the  test 
was  159  days.  The  test  was  carried  out  to  failure. 
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A  simple  creep  machine  was  designed  and  constructed  for 
the  test,  which  was  one  of  a  series  of  tests.  A  hydraulic 
system  with  a  bladder-type  accumulator  to  maintain  the 
necessary  load  constant  was  chosen.  The  loading  frame 
consisted  of  two  305  mm  x  305  mm  x  38  mm  steel  plates  spaced 
305  mm  apart  by  four  19  mm  high  tensile  steel  bolts  giving  a 
load  capacity  of  223  kN.  The  hydraulic  ram  is  ENERPAC  RC  256 
25  ton  cylinder  with  152  mm  stroke.  An  ENERPAC  P-39  single 

speed  hand  pump  drove  the  ram  and  pumped  up  the  accumulator. 

3 

The  accumulator  made  by  American  Bosch  is  328  cm  in  volume 
and  limits  the  system  pressure  to  52  MPa.  A  34.5  MPa  Marcsh 
pressure  gauge  of  0.25  percent  accuracy,  and  a  0.00254  mm 
(0.0001  inch)  dial  indicator  completed  the  required 
instrumentation.  The  creep  test  was  carried  out  at  room 
temperature  and  uncontrolled  humidity  and  the  sample  was 
loaded  perpendicular  to  the  bedding  planes.  The  laboratory 
set  up  is  outlined  in  Figure  3.8. 

The  data  from  the  creep  test  was  analyzed  as  outlined 
in  Section  3.3.  A  power  law  was  fitted  to  the  data  and  the 
data  points  and  the  best  fit  straight  lines  that  were  fitted 
to  the  decelerating  and  accelerating  creep  portions  were 
plotted  in  the  log  (strain  rate)  versus  log  (time)  plot.  As 
seen  in  Figure  3.9,  the  data  seem  to  scatter  but  they  fit 
well  to  the  power  law.  Figure  3.10,  a  plot  of  strain  rate 
versus  time,  shows  that  the  strain  rate  decreases  from  the 
beginning  to  a  certain  time  and  then  increases  towards 
failure.  Figure  3.10  clearly  illustrated  the  two  stages  of 
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Figure  3.8  Outline  of  Jeremic's  laboratory  set  up 
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Axial  Stress  =  19.4  MPa 


Figure  3.9  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  Starkey5 
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creep,  decelerating  and  accelerating. 

The  straight  line  fit  of  the  decelerating  creep  rate 

portion  in  Figure  3.9  has  a  slope  of  -1.00  and  a  one  minute 

strain  rate  (strain  rate  at  one  minute)  of  3971 

micro-strains  per  minute.  The  accelerating  creep  rate 

portion  has  a  slope  of  +8.09  and  a  strain  rate  value  at  one 
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minute  of  2.51x10  micro-strain  per  minute. 

Table  3.2  presents  the  results  of  analyzing  data  from 
the  Star-Key  Coal  Mine.  DW  stands  for  the  Durbin-Watson 
statistic  (Durbin  and  Watson,  1951)  which  is  a  test  for 
serial  correlation  in  the  residuals  of  the  fit.  If  the  creep 
law  is  a  reasonable  fit  to  the  data,  then  the  residuals  will 
be  randomly  distributed.  If  the  values  of  the  residuals  show 
some  dependence  on  the  variables,  x,  y,  then  the  proposed 
law  is  not  a  satisfactory  fit  to  the  data,  because  there 
still  remains  in  the  residuals  a  systematic  variation  which 
the  creep  law  has  not  satisfied.  Durbin  and  Watson  (1951) 
tabulated  two  groups  of  critical  values  for  DW  against  n, 
the  number  of  observations,  at  three  different  confidence 
levels.  The  5  per  cent  confidence  level  is  used  in  this 
analysis.  The  two  groups  of  critical  values  for  DW  are  the 
upper,  dU,  and  the  lower,  dL,  bounds.  If  the  observed  DW  is 
less  than  dL,  it  suggests  that  positive  serial  correlation 
of  the  residuals  exists  in  the  sample.  If  the  upper  value, 
dU,  is  not  exceeded,  positive  correlation  of  the  residuals 
might  exist  in  the  observations. 
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Table  3.2 


Fit 

Dece 1 erat Ing 


Accelerat Ing 


Summary  of  the  Star-Key  data  analysis 


Data 

Range 

log  A 
Intercept 

B 

Slope 

DW 

Slope 

Sign  If Icance 

log  Inflexion 
Time  (min) 

R 

Total  Fit 

1-80 

♦3.7564 

- 1 .0406 

1 .476 

73.945 

1-79 

+3.8045 

-1 .0526 

1 .459 

72.937 

1-78 

+3.9474 

-1 .0880 

1 .445 

76 . 196 

1-77 

+4.0498 

-1 . 1135 

1 .237 

73.891 

1-76 

+3.7616 

-1.0414 

1 .500 

77 . 274 

1-75 

+3.5991 

-1 .0007 

1 .677 

76.231 

1-74 

+3.5069 

-0.9776 

1 .779 

76 . 906 

1-73 

+3.5027 

-0.9765 

1 .759 

72 . 359 

1-72 

+3.3460 

-0.9373 

1.736 

63.208 

1-71 

+3.2557 

-0.9144 

1.721 

78.852 

60-90 

-41 .7652 

+7 . 7399 

1.617 

90.227 

5.190 

48 . 796 

61-90 

-44.8557 

+8.3276 

2.630 

112.218 

5 . 194 

49.215 

62-90 

-49.2026 

+9 . 1534 

1  .609 

197.507 

5.200 

49.609 

63-96 

-44.0275 

+8 . 1738 

1.681 

57 . 144 

5. 191 

49.039 

64-96 

-43.5989 

+8.0928 

1 .690 

52  815 

5 . 190 

48.973 

DW  =  Durbin  Watson  statistics 
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Axial  Stress  =  19.4  MPa 


Figure  3.10  Plot  of  strain  rate  (micro-strain/min )  vs  time 
(min)  Test  Starkey5 
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Another  measure  of  goodness  of  fit  is  the  Test  of  Slope 
Significance.  Cruden  (1971a)  outlined  the  calculation  of  R1  , 
the  ratio  of  the  estimated  variance  of  a  simple  linear 
regression  to  the  variance  about  the  mean  value  of  the 
dependent  variable.  The  statistic,  R1,  can  be  referred  to 
F-tables  with  one  and  (n-2)  degrees  of  freedom  with  1  per 
cent  confidence  level  suggested  by  Cruden  (1971a)  because  of 
the  large  amount  of  data. 

The  upper  bound  of  the  Durbin-Watson  statistic  is 
calculated  in  the  Subroutine  BFIT  in  package  CPACK.  The 
regression  analysis  will  carry  on  until  the  observed  DW 
exceeded  the  upper  bound  value  and  R1  exceeded  10. 

In  Table  3.2,  several  fits  of  the  decelerating  and 
accelerating  regression  lines  are  summarized.  For  the 
decelerating  fit,  the  regression  lines  fitted  to  data  range 
from  1  to  75  or  less  satisfied  the  two  criteria  outlined 
above.  The  rest  of  the  decelerating  fits  had  the  observed  DW 
too  low  to  exceed  the  upper  value,  dU.  The  regression  line 
fitted  to  data  range  from  1  to  75  was  chosen  as  the 
decelerating  fit.  For  the  accelerating  fits,  all  five 
regression  lines  summarized  in  Table  3.2  satisfied  the  two 
criteria.  It  appeared  that  the  regression  lines  with  data 
range  of  61  to  96  and  62  to  96  were  the  optimum  fits. 
However,  the  regression  line  with  data  range  of  64  to  96  was 
chosen  as  the  accelerating  regression  line  plotted  in  Figure 
3.9  because  it  was  the  smallest  range  of  data  used  to 
calculate  the  goodness  of  overall  fit,  R,  by  combining  the 
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decelerating  and  accelerating  fits  in  Subroutine  CFIT.  The 
inflexion  time  in  Table  3.2  is  the  time  when  the 
decelerating  regression  line  with  data  range  of  1  to  75 
intersects  with  the  accelerating  regression  line.  The 
inflexion  time  was  calculated  by  equating  the  decelerating 
and  accelerating  power  laws. 

The  scatter  in  the  Star-Key  Data  can  be  explained  by 
the  change  of  ambient  temperature  and  humidity.  A  very 
slight  change  in  temperature  and  humidity  can  affect  the 
deformation  of  a  sample  under  stress.  The  temperature  and 
humidity  dependent  deformation  are  being  recorded  together 
with  the  real  creep  deformation  and  this  causes  the  scatter 
of  the  data.  The  analysis  of  the  Star-Key  Data  is  included 
in  the  User's  Manual  in  Appendix  A  as  an  example. 

Finally,  base  on  the  analysis  of  the  experimental  data, 
the  form  of  the  creep  of  the  Star-Key  Coal  is  as  follow: 

-1.00  -44  +8.09 

e  =  397 1 t  +  (2.51x10  )t  (3.7) 

where  e  is  the  strain  rate  in  micro-strains  per  minute 
t  is  the  elapsed  time  in  minutes 


4.  Presentation  of  Test  Results 


4.1  Typical  Results 

As  mentioned  in  Section  2.3,  there  were  two  different 
sources  of  coal  samples  used  in  this  study.  The  first  source 
was  samples  obtained  from  the  sampling  site  and  the  second 
source  was  samples  not  used  by  da  Fontoura  (1980).  In 
Section  2.5,  the  three  different  types  of  tests  carried  out 
in  this  experimental  program  were  described  in  detail.  The 
first  type  was  testing  with  both  confining  and  back 
pressure,  referred  as  the  C-tests;  the  second  type  was  the 
uniaxial  compresison  test,  referred  as  the  U-tests;  and 
lastly,  testing  on  saturated  samples,  referred  as  the 
S-tests . 

Because  of  the  similarity  in  test  results,  not  all  the 
results  from  all  tests  are  presented  here.  Only  a  few  tests 
from  each  group  of  tests  are  presented  as  typical  results. 
The  rest  of  the  results  are  included  in  Appendix  B. 

Figure  4.1  to  Figure  4.7  are  double  logarithm  plots  of 
strain  rate  versus  time  of  the  results  of  tests  done  on 

samples  from  the  first  source.  Figure  4.1  and  Figure  4.2  are 

results  of  the  C-tests,  Figure  4.3  to  Figure  4.5  are  results 
of  S-tests  and  Figure  4.6  and  Figure  4.7  are  of  the  U-tests. 
Table  4.1  is  a  summary  of  the  decelerating  fit  of  all  the 
tests  carried  out  with  samples  from  the  first  source. 

Figure  4.8  to  Figure  4.16  are  double  logarithm  plots  of 

strain  rate  versus  time  of  the  results  of  tests  done  on 
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Table  4.1  Summary  of  decelerating  fit  to  data  from  tests 
using  fresh  samples 


Test 

(MPa) 

Stress 

log  A 
Intercept 

B 

Slope 

ConA 

ConB 

DW 

Slope 

S ign i f icance 

(min) 

tm 

C9-9*8 

3.99 

1 .903 

-1 .264 

0.  149 

0.053 

1 .734 

572.615 

366 

C9-9#9 

4.25 

1 .592 

-1.124 

0.  129 

0.  122 

1.811 

85.460 

10 

U9-9#9 

6 . 38 

1 .504 

-1 .059 

0.097 

0.034 

1 .758 

962 .546 

383 

U9-9#1 1 

6 . 38 

1 .607 

-1 .373 

0.059 

0.030 

2.375 

2035.032 

52 

U9-9/M7 

3 . 99 

1 .840 

-1.131 

0.114 

0.043 

2.417 

688.650 

240 

U9-9#21 

6 . 38 

1.817 

-1 .447 

0.131 

0.066 

1.671 

481 .304 

52 

U10-9*2 

7 . 18 

1 .844 

-1 .313 

0.  154 

0.053 

1 .914 

615.810 

478 

S9-9#7 

4 . 78 

1 .887 

-1 .077 

0.064 

0.029 

1 .649 

1360.900 

77 

S9-9#21 

4 . 78 

2. 122 

-1.175 

0.093 

0.059 

1 .585 

396.513 

22 

S9-9#9 

4.25 

1 . 705 

- 1 . 065 

0.070 

0.023 

1 . 704 

2207.271 

857 

S9-9#  1  1 

4 . 25 

1 .850 

-  1 . 122 

0.114 

0.038 

1 .696 

874 . 360 

620 

S9-9/M5 

3 . 99 

2.092 

-  1 .401 

0.098 

0.047 

1 .677 

889 . 559 

59 

S9-9#17 

3.99 

1 .977 

-1 .234 

0.088 

0.029 

1 .707 

1787.414 

597 

S 10-9^2 

3.99 

2 . 167 

-1 .344 

0.072 

0.025 

1.691 

2967.781 

461 

Con  =  confidence  limits  on  following  parameter 
DW  =  Durbin  Watson  statistics 


tm  *  Man  of  the  logarithm  of  time 
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Axial  Stress  =  3.99  MPa 


Figure  4.1  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  C9-9#8 


LOG  STRAIN  RATE 
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Axial  Stress  =  4.25  MPa 


Figure  4.2  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  C9-9#9 
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Axial  Stress  =  4.25  MPa 


Figure  4.3  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  S9-9#9 
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Axial  Stress  =  4.25  MPa 

Figure  4.4  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  S9-9#11 


LOC  STRAIN  RATE 


54 


LOG  TTTE 


Axial  Stress  =  4.78  MPa 


Figure  4.5  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  S9-9#21 
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Axial  Stress  =  6.38  MPa 

Figure  4.6  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  U9-9#9 
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Axial  Stress  =  3.99  MPa 


Figure  4.7  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  U9-9#17 
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samples  from  the  second  source.  Figure  4.8  and  Figure  4.9 
are  results  from  the  U-tests.  Figure  4.10  to  Figure  4.16  are 
plots  from  the  results  of  the  C-tests.  There  were  no  S-tests 
carried  out  with  samples  from  the  second  source.  Table  4.2 
presents  a  summary  of  the  decelerating  fit  of  all  the  tests 
carried  out  with  samples  from  the  second  source. 

In  the  decelerating  fit  of  a  regression  line  to 
experimental  data  presented  in  Figure  4.1  to  Figure  4.16, 
one  can  observe  that  as  time  increases  the  data  scattered 
more  about  the  regression  line.  This  is  due  to  the  fact  that 
the  observations  were  affected  by  the  change  of  the  ambient 
temperature.  At  the  beginning  of  the  test,  the  strain  rate 
is  high  and  therefore  not  significantly  affected  by  the 
temperature  effect.  With  the  strain  rate  decreasing,  it 
comes  to  a  point  where  the  temperature  effects  becomes 
significant  and  this  explains  the  scattering  of  data  as  time 
increases . 

All  except  one  of  the  double  logarithm  plots  presented 
here  are  decelerating  fits  of  the  test  results.  The  result 
from  Test  U9~9#9  is  the  only  test  that  showed  accelerating 
creep,  as  illustrated  in  Figure  4.6.  Table  4.3  summarized 
the  parameters  of  the  accelerating  fit.  Test  U9-9#9  was 
carried  out  with  a  relatively  high  stress  and  long  duration. 
There  were  other  tests  with  approximately  the  same  axial 
stresses  but  shorter  durations,  therefore  there  were  no 
accelerating  creep.  The  axial  stress  and  duration  of  all 
tests  are  summarized  in  Table  2.2. 
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Table  4.2  Summary  of  decelerating  fit  to  data  from  tests 
using  da  Fontoura’s  samples 


Test 

(MPa) 

Stress 

log  A 
Intercept 

B 

Slope 

ConA 

ConB 

DW 

Slope 

Sign if icance 

(min) 

tm 

UT6C 

2.92 

1 .715 

-1 .089 

0.040 

0.016 

1.701 

4444 . 102 

152 

UT6F 

1 . 86 

1 .738 

-1.041 

0.047 

0.020 

1.714 

2798.874 

124 

UT6F2 

2.39 

1 .976 

-1 .082 

0.084 

0.026 

2.296 

1773.229 

1 140 

UT6D 

2 . 39 

1 .578 

-0.913 

0.062 

0.022 

2.053 

1652 . 1 10 

338 

UT6D2 

2 . 39 

1.7  17 

-0.940 

0.077 

0.027 

2.414 

1218.440 

444 

UT6D3 

2.92 

1 . 865 

-1 .082 

0.  180 

0.051 

1 .732 

459.033 

2891 

CT6A1 

1  .  86 

1.818 

-1.125 

0.069 

0.024 

1 .703 

2119. 307 

372 

CT6A2 

2.66 

1 .269 

-0.958 

0.115 

0.037 

2.064 

657 . 122 

748 

CT6A3 

2.66 

1 . 754 

-1.186 

0.  127 

0.042 

1  .691 

783 . 276 

57  1 

CT6A4 

2.66 

1 . 559 

-0.919 

0. 195 

0. 154 

1  .608 

35.451 

14 

CT6A5 

3 . 72 

1 .646 

-0.955 

0.065 

0.02  1 

1.917 

2 1 14.035 

843 

CT6A6 

3.72 

1.619 

-1 .008 

0.049 

0.018 

1.691 

3069 . 792 

246 

CT6A7 

3.72 

1 .779 

-  1  .04  1 

0.099 

0.035 

2 . 294 

880.733 

379 

CT6A8 

4.78 

1.717 

-0.941 

0.054 

0.019 

1  .921 

2436 . 323 

427 

CT6A9 

5 . 32 

1.811 

-1.043 

0.054 

0.021 

1  .696 

2490.586 

238 

Con  =  confidence  limits  on  following  parameter 
DW  =  Durbin  Watson  statistics 
tm  *  mean  of  the  logarithm  of  time 
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Table  4.3  Summary  of  accelerating  fit 


U9-9#9 

Acce 1 erat Ing  Fit 

Stress  (MPa)  6.38 

Intercept,  Log  A  -19.333 

Slope,  B  4.272 

ConA  3.492 

ConB  0.888 

DW  6.401 

Slope  Significance  23.180 

tm  (min)  8590 

Con  *  confidence  limits  on  following  parameter 
DW  =  Durbin  Watson  Statistics 
tm  =  mean  of  the  logarithm  of  time 
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Axial  Stress  =  2.92  MPa 


Figure  4.8  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  UT6C 
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LX  TIHE 


Axial  Stress  =  1.86  MPa 


Figure  4.9  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  UT6F 
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Axial  Stress  =  2.66  MPa 


Figure  4.10  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A2 
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Axial  Stress  =  2.66  MPa 

Figure  4.11  Logarithm  plot  of  strain  rate  (micro~stra in/min ) 
vs  time  (min)  Test  CT6A4 
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Axial  Stress  =  3.72  MPa 

Figure  4.12  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A5 
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Axial  Stress  =  3.72  MPa 

Figure  4.13  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A6 
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Axial  Stress  =  3.72  MPa 

Figure  4.14  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A7 
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Axial  Stress  =  4.78  MPa 


Figure  4.15  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A8 
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Axial  Stress  =  5.32  MPa 

Figure  4.16  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A9 
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4.2  Comparison  with  da  Fontoura's  Test  Results 

da  Fontoura  (1980)  reported  nine  single  stage  creep 
tests  on  Wabamun  Coal.  The  method  of  data  analysis  was 
similar  to  the  one  mentioned  in  Section  3.3,  except  that 
there  was  no  computer  program  employed  to  evaluate 
accelerating  creep.  Table  4.4  summarizes  da  Fontoura's 
results  of  nine  single  stage  creep  tests. 

All  of  da  Fontoura's  tests  summarized  in  Table  4.4  were 
confined  tests,  with  confining  pressure  ranging  from  208  kPa 
to  553  kPa  (da  Fontoura,  1980).  The  effect  of  confining 
pressure  is  to  decrease  the  strain  per  cracking  event  and  so 
decreases  strain  rate  and  postpones  the  onset  of 
accelerating  creep  (Lama,  1978,  pp. 251-253).  The  test 
results  summarized  in  Table  4.1  are  from  tests  carried  out 
with  samples  which  are  of  a  different  orientation  to  da 
Fontoura's  samples.  Therefore,  the  results  in  Table  4.1  and 
Table  4.4  can  be  compared  only  qualitatively. 

Results  summarized  in  Table  4.2  were  from  tests 
performed  on  samples  from  the  second  source,  samples  not 
used  by  da  Fontoura  (1980).  Comparing  the  results  of  the 
confined  tests  in  Table  4.2  to  results  in  Table  4.4,  the 
results  in  Table  4.2  had  higher  one  minute  strain  rate 
(strain  rate  at  one  minute)  values  and  slightly  steeper 
slopes,  da  Fontoura  (1980)  showed  that  the  one  minute  strain 
value  of  the  regression  analysis  was  stress  level  dependent 
but  there  was  no  indication  of  any  relationship  between  the 
slope  value  and  the  stress  level.  One  can  argue  that  samples 
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Table  4.4  Summary  of  regression  analysis  on  single  stage 
creep  tests  reported  by  da  Fontoura 


Test 

(MPa) 

Dev.  Stress 

log  A 
Intercept 

B 

S  1  ope 

CT  1 

0.70 

1 .439 

-0.896 

CT2 

0.50 

1 .380 

-0.856 

CT3 

1 .60 

1 .274 

-0.819 

CT  4 

2.20 

1  .473 

-0.882 

CT6 

2 . 18 

1 .479 

-1 .040 

CT7/St 1 

2.88 

1.417 

-0.931 

CT7/St2 

6.00 

2.322 

-0.810 

CT8 

3.57 

1 .428 

-0.919 

CT9 

4.02 

1 .566 

-0.994 

. 
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not  used  by  da  Fontoura  were  desiccated,  even  though  they 
were  stored  in  the  moisture  room,  however  for  a  considerably 
long  period  of  time  (about  three  years).  There  were  more 
cracks  observed  in  these  samples  than  in  the  fresh  samples. 
As  the  samples  were  weakened  by  desiccation,  the  stress 
level  would  increase  because  of  the  reduced  strength.  This 
might  be  the  reason  for  the  high  one  minute  strain  values 
reported  in  Table  4.2  as  compared  to  da  Fontoura' s  results 
in  Table  4.4.  The  higher  values  of  B  shown  in  Table  4.2  are 
probably  due  to  the  closure  of  cracks  as  load  is  applied. 

4.3  Results  from  Additional  Tests 

Two  more  tests  were  carried  out  using  Jeremic's 
(Personal  Communication)  laboratory  apparatus,  which  were 
described  in  detail  in  Section  3.4.  The  tests  were  labelled 
J9-9#9  and  J9  —  9# 11,  J  stands  for  the  type  of  tests  using 
Jeremic's  apparatus. 

There  is  a  slight  difference  between  the  laboratory  set 
up  and  Jeremic's  original  apparatus.  A  LVDT  is  used  instead 
of  a  dial  gauge  to  monitor  the  deformation  of  the  sample. 

The  LVDT  is  a  24-DCDT-250  LVDT  manufactured  by  Hewlett 
Packard,  with  a  displacement  range  of  ±6.35  millimetres.  The 
LVDT  was  clamped  on  to  the  hydraulic  ram  and  the  advance  of 
the  ram  was  recorded  as  the  axial  deformation  of  the  sample. 
A  power  supply  unit  capable  of  providing  input  voltage  of  24 
volts  was  used  to  feed  the  LVDT.  The  displacement  was 
monitored  by  the  Hewlett  Packard  3054  data  logger  and  all 
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the  observations  were  stored  on  a  disc  with  the  use  of  a 
HP82901M  Flexible  Disc  Drive  Unit. 

The  tests  were  uniaxial  compression  test  with  a 
constant  load  applied  by  a  hydraulic  ram  driven  by  a  hand 
pump.  For  Test  J9-9#9,  the  axial  stress  was  10.65  MPa.  The 
test  lasted  for  approximately  twelve  days  before  rupture 
occurred.  The  axial  stress  for  Test  J9— 9# 1 1  was  9.74  MPa  and 
the  test  lasted  for  about  fourteen  days  until  rupture.  In 
both  tests,  there  were  slight  bulging  in  the  sample  and  some 
spalling  from  the  sample  as  the  test  went  on.  This  will 
result  in  a  reduction  of  the  cross  section  of  the  sample 
taking  the  load  and  a  consequent  increase  in  the  stress 
acting  on  the  load  bearing  portion  of  the  sample.  It  is 
possible  then  that  these  tests  are  not  true  creep  tests, 
that  is  they  were  not  carried  out  under  constant  stress. 

The  double  logarithm  plots  of  strain  rate  versus  time 
of  the  results  from  Test  J9-9#9  and  Test  J9— 9# 1 1  are 
presented  in  Figure  4.17  and  Figure  4.18  respectively.  The 
regression  analysis  for  both  tests  are  summarized  in  Table 
4.5. 

For  the  results  from  Test  J9— 9#1 1 ,  fifteen  observations 
had  to  be  truncated  from  the  end  of  the  input  file  before  an 
accelerating  creep  law  could  be  fitted  to  the  data. 
Observations  were  taken  once  every  two  hours  at  the  end  of 
the  test,  fifteen  observations  would  be  equivalent  to  thirty 
hours  in  actual  time.  The  reason  for  the  truncation  of  the 
data  is  that  there  was  too  much  scatter  towards  the  end  of 
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Table  4.5  Summary  of  fits  to  Test  J9-9#9  and  Test  J9-9#11 


J9-9#9  J9-9*11 


Decelerat ing 

Acce 1 er at i ng 

Decelerat ing 

Accelerat ing 

Stress  (MPa) 

10.65 

1065 

9 . 74 

9 . 74 

Intercept,  log  A 

2 . 166 

-12.841 

1 .722 

-93 . 139 

Slope,  B 

-0.770 

3.053 

-0.492 

21 .737 

ConA 

0.110 

1  .628 

0.030 

26.448 

ConB 

0.092 

0.405 

0.01  1 

6 . 188 

DW 

1 .662 

1 .705 

1 .763 

1 .942 

Slope  Significance 

70.098 

53.783 

2102.465 

10.552 

tm  (min) 

1  1 

10520 

364 

18793 

Total  Fit,  R 

388.931 

380.018 

Con  =  confidence  limit  on  following  parameter 
DW  =  Durbin  Watson  statistics 
tm  =  mean  of  the  logarithm  of  time 
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Axial  Stress  =  10.65  MPa 


Figure  4.17  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  J9-9#9 
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Axial  Stress  =  9.74  MPa 


Figure  4.18  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  J9-9#11 
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the  test,  which  creates  a  problem  for  the  least  squares 
calculation.  The  scatter  is  most  likely  due  to  the 
fluctuation  in  deformation  caused  by  the  falling  of  spalls 
of  material,  which  will  result  in  a  reduction  of  the  cross 
section  of  the  sample  and  a  consequent  increase  in  the 
stress  and  deformation  rate.  The  spalling  observed  in  Test 
J9-9#11  was  more  severe  than  in  Test  J9-9#9.  This  would 
indicate  that  the  stress  and  deformation  rate  increase  in 
Test  J9-9#11  was  higher  and  the  consequent  fluctuation  in 
deformation  was  higher  too. 

The  variability  of  the  material  is  also  an  important 
factor  to  the  deformation  properties.  Samples  from  the  same 
sampling  site  may  have  very  different  deformation  properties 
(Kaiser  and  Maloney,  1982),  and  this  is  reflected  in  the 
experimental  observations.  The  variability  of  the  material 
is  shown  by  the  results  from  Test  J9-9#9,  where  no  data 
truncation  is  required  for  a  good  fit. 


5.  Final  Remarks 


5.1  Discussion  of  Test  Results 

The  test  results  from  each  laboratory  test  were  fitted 
to  a  decelerating  power  law  and  the  decelerating  fits  to  all 
the  tests  were  summarized  in  Table  4.1  and  Table  4.2.  The 
results  in  Table  4.1  were  from  tests  carried  out  with  fresh 
samples  from  the  sampling  site,  while  Table  4.2  contained 
results  from  tests  performed  on  samples  not  used  by  da 
Fontoura  (1980). 

da  Fontoura' s  samples  were  drilled  from  a  block  sample 
with  their  long  axis  parallel  to  their  bedding  planes  and  at 
an  angle  of  30°  with  the  major  cleat.  This  orientation  would 
give  the  sample  the  minimum  axial  strength  due  to  the 
presence  of  the  major  cleat  at  30°  to  the  long  axis  (Hoek 
and  Brown,  1980,  pp.  157-165),  and  it  would  also  cause 
pre-existing  cracks  that  were  parallel  to  the  cleavage  to 
rupture  in  shear  when  loaded  axially  in  compression.  For 
samples  that  were  drilled  with  their  long  axes  perpendicular 
to  their  bedding  planes,  cracks  that  were  parallel  to  the 
bedding  would  be  closed  but  would  not  rupture  in  shear  as 
load  is  applied. 

By  comparing  the  C-test  results  in  Table  4.1  and  Table 
4.2,  one  can  calculate  the  strain  rates  by  assigning  a  log  t 
value  to  the  power  laws  and  find  out  that  the  C-tests  in 
Table  4.2  are  creeping  at  a  faster  rate  than  those  in  Table 

4.1  as  time  increases.  For  example,  compare  test  C9-9#8  from 
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Table  4.1  to  test  CT6A6  from  Table  4.2.  When  (log  t)  is  1.0, 
(log  e )  for  test  C9-9#8  and  CT6A6  are  0.639  and  0.611 
respectively.  When  (log  t)  is  2.0,  (log  e)  for  test  C9-9#8 
is  -0.625  as  compared  to  -0.397  for  test  CT6A6.  This  showed 
that  samples  drilled  with  their  axes  at  30°  to  the  joint  are 
creeping  more  rapidly  than  those  samples  drilled  with  their 
axes  perpendicular  to  the  bedding  planes.  This  is  the  kind 
of  result  that  one  would  expect  as  a  result  of  difference  in 
orientation  in  the  two  types  of  samples. 

The  slope  values  in  Table  4.1  were  higher  than  those 
reported  in  Table  4.2.  This  could  probably  be  explained  by 
the  difference  in  orientation  of  the  samples.  The  slope 
value  is  a  strain-hardening  parameter  measuring  the  rate  of 
the  decrease  of  strain  rate  with  time.  Cruden  (1970) 
suggested  that  for  slope  values  that  are  less  than  -1,  the 
creep  strains  will  approach  a  finite  value.  For  slope  values 
that  are  greater  than  -1,  the  creep  strain  will  increase  to 
a  critical  creep  strain  where  accelerating  creep  begins. 
There  was  only  one  test.  Test  U9-9#9,  showed  accelerating 
creep.  This  test  has  a  slope  of  -1.059  for  the  decelerating 
fit,  which  is  the  least  value  of  slope  in  Table  4.1.  Taking 
the  confidence  limit  of  the  slope  value  into  account.  Test 
U9-9#9  could  have  a  slope  value  greater  than  -1. 

Test  U9-9#9  was  the  only  test  in  which  an  accelerating 
creep  regression  line  could  be  fitted  to  the  data.  The 
accelerating  creep  shown  was  not  as  prolonged  as  the 
Star-Key  Coal  test  because  the  test  was  not  carried  out  to 
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rupture . 

In  the  C-tests,  the  application  of  confining  pressure 
restricted  the  lateral  displacement  of  the  samples  and  they 
became  stronger  (Jaeger  and  Cook,  1969,  pp.  86-88). 

Confining  pressure  also  decreases  strain  rate  and  strain  per 
cracking  event,  which  in  turn  postpones  the  onset  of  the 
accelerating  creep. 

Hoek  (1965)  showed  that  rock  strengths  are  moisture 
sensitive.  The  strength  of  rock  is  reduced  when  'wet'  and 
the  influence  of  moisture  is  more  pronounced  in  the  case  of 
materials  such  as  coal  (Hoek,  1965,  pp.  118-122).  In  Table 
4.1,  results  from  the  S-tests  did  not  show  any  distinct 
difference  from  other  type  of  tests.  This  may  be  due  to  the 
fact  that  samples  for  the  S-tests  were  submerged  in  water 
for  only  24  hours  prior  to  testing,  while  in  Hoek's  report, 
'wet'  referred  to  samples  submerged  in  water  for  120  days. 

Cruden  (1974)  developed  a  static  fatigue  law  which 
relates  the  strength  of  the  material  to  the  time  that  it  had 
been  under  stress.  The  form  of  the  law  depends  on  the 
relationship  between  stress  and  creep  rate,  i.e.,  the  form 
of  the  stress  dependence  of  the  creep  rate.  In  an  attempt  to 
find  the  form  of  the  stress  dependence  of  the  creep  rate  of 
coal,  tests  were  grouped  by  the  source  of  samples  and  then 
grouped  again  by  the  type  of  tests.  As  the  slope  values  of 
the  decelerating  fits  from  most  of  the  test  results  were 
less  than  -1  and  a  wide  range  of  stresses  was  not  used  in 
each  group  of  tests,  the  form  of  stress  dependence  of  the 
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creep  rate  could  not  be  found. 

The  same  static  fatigue  law  was  applied  to  test  results 
reported  by  da  Fontoura  (1980).  The  form  of  stress 
dependence  of  the  creep  rate  could  not  be  found  either 
because  the  points  on  the  double  logarithm  plot  of  the  ratio 
of  the  strain  rates  versus  the  ratio  of  the  stresses  were 
too  scattered  for  a  significant  linear  regression. 

The  test  results  presented  in  Table  4.1  and  Table  4.2 
all  showed  a  slope  value  of  the  decelerating  fit  of  less 
than  -1.  This  prevented  the  calculation  of  the  critical 
strain  for  the  onset  of  the  accelerating  creep.  Tests 
carried  out  with  Jeremic’s  apparatus,  Test  J9-9#9  and  Test 
J9-9#11  were  the  only  tests  that  had  an  accelerating  fit  and 
the  slope  of  the  decelerating  fit  greater  than  -1.  By 
integrating  the  decelerating  power  law  and  taking  the 
intersects  of  the  decelerating  and  accelerating  fitted  lines 
as  the  time  when  accelerating  creep  started.  The  critical 
strains  for  Test  J9-9#9  and  Test  J9— 9# 1 1  were  calculated  as 
1.5251  and  1.5276  percent  strains,  respectively. 

In  this  research  program,  the  variability  of  the  coal 
represents  a  major  source  of  problems.  From  the  results  of 
the  sample  characterization  presented  in  Table  2.1,  one  can 
calculate  the  mean  and  the  standard  deviation  of  the 
compressional  elastic  wave  velocities,  Vp.  The  standard 
error  of  the  mean  can  be  calculated  by  dividing  the  standard 
deviation  by  the  mean.  The  standard  error  of  the  mean  of  the 
compressional  elastic  wave  velocities  of  the  fresh  samples 
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is  calculated  as  0.1378  and  for  da  Fontoura' s  samples  is 
0.0505.  Cruden  (1969)  listed  values  of  the  compressional 
elastic  wave  velocities  of  samples  of  Carrara  Marble  and 
Pennant  Sandstone.  The  standard  error  of  the  mean  of  the 
compressional  elastic  wave  velocities  of  Carrara  Marble  is 
0.0384  and  for  Pennant  Sandstone  is  0.0075.  By  comparing 
these  figures,  one  can  conclude  that  coal  is  a  more  variable 
material  than  Carrara  Marble  and  Pennant  Sandstone. 

Kaiser  and  Maloney  (1982)  investigated  the  deformation 
properties  of  a  Sub-bituminous  coal  mass  from  the  same  site 
by  conducting  a  series  of  compression  tests  on  large  block 
samples  of  coal.  The  variability  of  the  Young's  Modulus,  E, 
of  the  coal  mass  was  between  0.95  GPa  to  3.25  GPa  for  one 
sample  and  between  0.45  GPa  to  2.90  GPa  for  the  other 
sample.  They  concluded  that  coal  is  a  highly  heterogeneous 
material  and  the  bulk,  K,  and  the  Young's  Modulus,  E,  can 
vary  widely  within  a  small  area.  Therefore,  one  can 
reasonably  argue  that  the  lack  of  obvious  stress  dependence 
of  the  creep  rate  in  the  test  results  is  due  to  the 
variability  of  the  coal  samples. 

da  Fontoura  (1980)  suggested  that  for  the  decelerating 
power  law,  the  one  minute  strain  rate  ,  A,  is  stress 
dependent  but  the  slope,  B,  is  essentially  independent  of 
the  stress  applied.  Using  the  test  results  of  the  single 
stage  creep  tests  reported  by  da  Fontoura  (1980)  in  Table 
4.4,  the  plots  of  (log  A)  versus  (stress)  and  B  versus 
stress  were  presented  in  Figure  5.1  and  Figure  5.2, 
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respectively. 

In  da  Fontoura’s  results,  the  result  from  creep  test 
CT7/st2  was  excluded  from  Figure  5.1  and  Figure  5.2  because 
it  was  the  second  stage  of  a  multiple  stage  creep  test.  In 
the  legend  for  both  Figure  5.1  and  Figure  5.2,  the  first 
letter  represents  the  type  of  test,  the  numbers  represent 
the  angle  between  the  sample  axis  and  the  major  joint  or  the 
bedding  planes.  The  last  two  letters  in  the  legend  stand  for 
the  initials  of  the  person  who  reported  the  tests.  The 
infilled  symbols  are  tests  with  sample  axes  perpendicular  to 
bedding  planes.  The  crosses  are  tests  carried  out  using 
Jeremic's  apparatus. 

In  Figure  5.1,  the  one  minute  strain  rate.  A,  has  a 
unit  of  micro-strains  per  minute.  There  were  no  obvious 
stress  dependence  of  A  after  the  exclusion  of  test  CT7/st2 
from  da  Fontoura's  results.  Results  from  laboratory  tests 
reported  in  Section  4.1  were  also  plotted  in  Figure  5.1, 
they  seemed  to  support  the  same  statement. 

The  plot  of  slope,  B,  versus  stress  is  presented  in 
Figure  5.2.  For  each  group  of  tests  represented  by  the 
legend,  there  seemed  to  be  no  direct  relationship  between 
the  slope  of  the  decelerating  power  law  and  the  stress 
applied. 

Therefore,  it  is  reasonable  to  suggest  that  for  a 
selection  of  coal  specimens  subject  to  single  stage  creep 
tests,  the  one  minute  strain  rate  and  the  slope  of  the 
fitted  decelerating  power  law  show  no  clear  dependence  on 
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Figure  5.1  Plot  of  (log  A)  vs  (stress) 
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Figure  5.2  Plot  of  (Slope,  B)  vs  (Stress) 
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the  stress  applied. 

For  tests  at  higher  stresses,  samples  with  parameter  B 
less  than  -1  are  creeping  under  stresses  below  their  long 
term  strength.  Examination  of  Figure  5.2  shows  that  the 
parameter  B  of  the  tests  represented  by  the  crosses  are 
greater  than  -1,  thus  they  are  creeping  under  stresses  above 
their  long  term  strength.  Test  U9-9#9  could  possibly  be 
creeping  under  a  stress  above  its  long  term  strength  because 
its  B  value  is  so  close  to  1  and  it  appears  to  have  shown 
accelerating  creep. 

5.2  Summary  and  Conclusion 

Coal  samples  were  obtained  from  the  Highvale  Mine, 
Alberta.  These  samples  together  with  those  not  used  by  da 
Fontoura  (1980)  were  used  for  the  creep  experiments.  The 
creep  rig  shown  in  Figure  2.3  was  used  and  all  the 
experimental  observations  were  conditioned,  processed  and 
analyzed  with  the  computer  package  CPACK. 

The  experimental  data  were  first  fitted  with  a 
decelerating  power  law  as  suggested  by  Cruden  (1971a),  and 
the  fitting  of  an  accelerating  power  law  was  also 
investigated.  The  test  results  from  the  Star-Key  Coal  test 
carried  out  by  Jeremic  (Personal  Communication)  showed  that 
it  is  possible  to  describe  the  accelerating  creep  with  an 
accelerating  power  law.  The  results  from  this  experimental 
program  were  summarized  in  Table  4.1  and  Table  4.2.  Most  of 
the  double  logarithm  plots  of  strain  rate  versus  time  were 
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shown  in  Figure  4.1  to  Figure  4.16  and  the  rest  of  the  plots 
were  put  together  in  Appendix  B. 

For  those  experiments  that  showed  accelerating  creep, 
the  use  of  an  accelerating  power  law  seemed  to  describe  the 
accelerating  creep  adequately.  Therefore  it  may  be  concluded 
that  the  entire  creep  curve  can  be  described  by  two  power 
laws  as  follow: 

B  D 

e  =  At  +  Ct  (3.4) 

The  attempt  to  find  the  stress  dependence  of  the  creep 
rate  of  a  model  coal  pillar  was  frustrated  by  the 
variability  of  the  material.  The  Young's  Modulus,  E,  of  coal 
can  vary  widely  within  a  relatively  small  area,  which 
indicates  that  samples  from  the  same  sampling  site  can  have 
much  different  deformation  properties.  This  suggests  that 
the  parameters  of  the  creep  curves  determined  in  the  tests 
cannot  be  reliably  be  scaled  up  for  design  purposes. 

5.3  Recommendations 

To  avoid  the  problem  of  material  variability,  one  may 
match  samples  by  the  value  of  compressional  elastic  wave 
velocity,  Vp,  from  a  much  larger  group  of  samples.  However, 
one  major  problem  is  that  sampling  in  coal  is  not  an  easy 
task . 

Another  alternative  is  to  perform  increment  tests. 
Cruden  (1983)  suggested  that  the  stress  dependence  of  creep 
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rate  can  be  determined  by  an  increment  test  on  a  single  rock 
specimen.  Such  a  test  eliminates  variation  caused  by 
variation  of  physical  properties  between  specimens. 

Further  research,  using  one  of  these  two  alternatives 
might  allow  the  determination  of  the  physical  parameters 
controlling  the  creep  of  coal.  However,  results  from  both 
alternatives  may  be  misleading.  The  heterogeneity  of  coal 
induces  stress  concentration  which  may  accelerate  creep,  and 
cause  failure  modes  which  may  not  appear  in  more  uniform 
spec imens . 
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PROGRAM:  CPACK 

This  program  analyzes  data  from  creep  experiments  to 
identify  both  the  decelerating  and  accelerating  components 
of  creep. 

CPACK  calls  the  Subroutine  CRED  to  reduce  the 
experimental  data  and  compute  engineering  strains. 

The  Subroutine  BFIT  is  then  called  to  convert  the 
strains  into  strain  rates,  and  to  obtain  a  best  fit  of  a 
power  law  to  the  initial  portion  of  the  strain  rate  data. 

The  best  fit  decelerating  strain  is  then  computed  in 
the  Subroutine  INTEG  by  integrating  the  power  law 
expressions  for  decelerating  strain  rate  obtained  in  BFIT. 
This  decelerating  strain  is  then  subtracted  from  the 
experimental  strain  to  leave  accelerating  strain. 

BFIT  is  called  again  to  fit  a  power  law  to  the 
accelerating  strains  computed  in  INTEG.  On  this  call  to 
BFIT,  a  fit  is  made  of  the  latter  portion  of  the  data. 

At  this  point  the  user  has  the  option  of  calling 
subroutine  CFIT,  which  generates  a  table  of  fit  parameters 
and  statistics  for  a  number  of  trial  fits  of  the 
accelerating  data.  This  allows  the  user  to  select  the 
optimum  overall  fit  of  a  power  law  expression  to  the  data. 

Subroutine  CPLOT  is  then  called  to  produce  two  plots  of 
the  data:  Time  vs  Strain  rate,  and  Log  Time  vs  Log  Strain 
rate . 

Capabilities 

>  The  maximum  number  of  records  of  data  is  1000. 
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>  The  program  can  accept  data  in  the  standard  form 
indicated  below,  or  in  user  specified  forms. 

>  The  maximum  number  of  LVDTs  is  two. 

>  User  responses  to  prompts  must  be  in  capital  letters,  or 
if  numbers,  must  be  terminated  with  a  comma. 

>  The  first  minutes  of  a  creep  experiment  may  show  strains 
that  reflect  changes  of  the  load  on  the  specimen.  These 
changes  may  cause  problems  in  the  program.  Judgement 
should  therefore  be  used  in  selecting  the  first  records 
to  be  processed. 

>  If  the  data  does  not  contain  accelerating  creep,  the 
program  will  likely  generate  error  messages  on  the 
second  call  to  BFIT.  If  this  is  anticipated,  the  user 
can  choose  to  plot  the  decelerating  creep  data  before 
the  program  attempts  to  isolate  accelerating  creep. 

>  The  graphs  that  subroutine  CPLOT  can  generate  are 
optional . 

INPUT 

UNIT  4  -  Disk  -  Input  Data  File 

The  standard  input  data  file  is  described  in 
Table  1,  although  other  formats  can  be  specified  by  the 
user  in  the  prompt  sequence. 

Record  length  is  109  bytes.  Twelve  variables  are 
read  in  each  record  with  the  format  (7X,  512,  IX,  15, 
213,  4G20 ) . 


. r  ^  1 


' 


This  file  is  output  from  an  HP  3054  Data  Logger. 


TABLE  1  -  INPUT  DATA  FILE  FORMAT 
VARIABLE  NO.  COLUMNS  SPECIFICATIONS 


1 

008  - 

009 

Integer , 

Time  -  Month 

2 

010  - 

011 

Integer , 

Time  -  Day 

3 

012  - 

013 

Integer , 

Time  -  Hour 

4 

014  - 

015 

Integer , 

Time  -  Minute 

5 

016  - 

017 

Integer , 

Time  -  Second 

9 

030  - 

049 

Real,  Load 

10 

050  - 

069 

Real,  Cell  Pressure 

1  1 

070  - 

089 

Real,  Defl.  Reading  LVDT 

#1 

12 

090  - 

109 

Real,  Defl.  Reading  LVDT 

n 

UNIT  5  -  Terminal-User  Responses  to  Interactive  Prompts 
The  user  is  prompted  for  the  sample  dimensions,  and 
the  input  data  file  format,  if  it  deviates  from  the 
standard  format.  If  the  plotting  option  is  chosen,  the 
user  is  asked  to  input  the  length,  scale,  and  origin  of 
each  axis,  for  the  Time  vs  Strain  Rate  graph.  The 
default  for  these  parameters  is  a  lOin  X  lOin  graph 
containing  the  entire  range  of  data. 


OUTPUT 


UNIT  6  -  Terminal  -  Interactive  Prompts 

UNIT  7  -  Printer 

The  maximum  record  length  is  84  bytes.  Sample 
dimensions  and  a  table  of  strain  and  time  are  output  on 
this  unit.  In  addition,  a  table  of  transformed  data,  and 
fit  parameters  are  printed  for  each  strain  gauge,  for 
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both  accelerating  and  decelerating  creep.  A  table  of  fit 
parameters  for  a  number  of  trial  fits  can  also  be 
output,  if  desired,  to  allow  the  user  to  select  the 
optimum  overall  fit. 

UNIT  9  -  Output  File  -  PDF 

PDF  is  the  Plot  Description  File.  This  is  the 
information  required  by  the  Calcomp  plotter  to  plot  the 
graphs . 

UNIT  98  -  Output  File  -  Summary 

Output  in  the  form  of  a  summary  for  each  plot  is 
contained  in  -  summary. 

STORAGE  REQUIREMENTS 


Code:  CPACK  Bytes:  57,  344 

Disc  Pages:  14 


SUBROUTINES 


CRED 

Reduces  experimental  creep  data,  and  computes 


engineering  strain 
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SOLUTION 

The  time  clock  readings  are  converted  to  elapsed  time 
in  hours. 

Engineering  microstrains,  ESTRN ( J ) ,  are  computed  from 
the  deflection  readings,  RDEF(J),  using  the  following 
algorithm: 


VDEF(J)  =  RDEF(J)  -  RDEF ( 1 ) 

ESTRN ( J )  =  ( VDEF ( J ) /XLEN )  X  1000000  X  FACT 

Where  XLEN  is  the  sample  length  and  FACT  is  the  LVDT 
calibration  factor. 

CRED  writes  the  sample  parameters,  initial  readings, 
and  a  table  of  time  and  engineering  strain. 

Calling  Sequence:  CALL  CRED  (DIAM,  XLEN,  NRR ,  El,  E2,  TT, 
XANS1,  IANS3,  I  SAMI ,  ISAM2,  EZPST,  ESTRN 1 ,  ESTRN2 ) 

DIAM  -  Sample  Diameter 
XLEN  -  Sample  Length 

NRR  -  Number  of  records  passed  on  to  BFIT 

XANS1  -  Interactive  response  re:  Is  data  in  Standard  Form? 
IANS3  -  Number  of  deflection  guages 
I  SAMI,  ISAM2  -  Sample  number 
EZPST  -  Elapsed  time 

ESTRN 1  -  Engineering  strain  for  LVDT 
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ESTRN2  -  Engineering  strain  for  LVDT  #2 

TT ,  El  and  E2  are  identical  to  the  above  three 
variables  respectively,  except  that  the  first  record  is 
deleted.  These  three  are  the  variables  passed  on  to  BFIT. 

BFIT 

Fits  a  power  law  to  the  experimental  data. 

SOLUTION 

The  strains  are  made  consistently,  increasing  by  using 
an  averaging  technique,  outlined  on  p.  113  (Cruden,  1971a) 
which  smooths  out  the  portions  of  the  data  where  the  strains 
decreases.  These  strains  are  then  converted  to  strain  rates. 
Logarithms  of  time  and  strain  rate  are  then  taken,  to  allow 
a  least  squares  linear  regressions  to  fit  a  power  law  to  the 
data . 

The  program  computes  the  Durbin  Watson  Statistic  and 
the  Test  of  Slope  Significance  (Cruden,  1971)  to  evaluate 
the  goodness  of  fit.  If  the  Durbin  Watson  Statistic  is  less 
than  the  tabulated  upper  limit,  or  if  the  Slope  Significance 
is  less  than  10,  the  last  record  is  deleted  and  the  fitting 
process  is  repeated  with  the  remaining  records.  This 
continues  until  the  statistical  limits  are  satisfied,  and  a 
fit  of  the  initial  portion  of  the  data  is  achieved.  This 
yields  a  power  law  expression  for  the  strain  rate  of 
decelerating  creep. 
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On  the  second  call  to  BFIT  the  data  contains 
accerelating  creep  strains,  output  from  subroutine  INTEG. 

The  program  proceeds  as  for  decelerating  creep  on  the 
first  call,  except  that  each  time  the  statistical  limits 
orient  satisifed,  it  deletes  the  earliest  record,  until  a 
satisfactory  fit  of  the  latter  portion  of  the  data  is 
achieved.  The  program  then  writes  the  remaining  data  records 
in  the  fit,  along  with  the  slope,  intercept,  and  other 
statistical  parameters,  including  the  confidence  limits  on 
the  slope  and  intercept. 

The  estimates  of  the  slope  and  intercept  follow  the 
statistical  t-distr ibut ions ,  so  to  obtain  the  percentile 
values  of  confidence  limits  of  these  estimates,  the  user 
should  consult  the  widely  available  t-distr ibution  tables, 
and  determine  the  percentile  value  that  corresponds  to  the 
confidence  limit  listed  in  the  output,  for  the  number  of 
degrees  of  freedom  associated  with  the  data.  A  fit  is  made 
for  each  strain  gauge. 

Calling  Sequence:  CALL  BFIT  (NRR,  El,  E2,  TT,  L,  N,  BBO, 

BB 1 ,  XANS1,  IANS3,  LT,  LE ,  W1 ,  EEM1 ,  NCI,  K,  00,  FF ,  DW, 

MMM ,  ANSM3 ,  NRF ) 

Variables  not  defined  previously: 

L  -  Indicates  type  of  data 

L  =  1  -  decelerating  creep 
L  =  2  -  accelerating  creep 
N  -  record  number 
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BBO  -  array  containing  intercepts  of  fitted  lines 
BB 1  -  array  containing  slopes  of  fitted  lines 

The  following  variables  are  either  passed  to,  or  received 
from  subroutine  CFIT: 

LT,  LE,  and  W1  -  arrays  that  contain  log  time,  log  strain 
rate,  and  the  record  weighting  respectively,  calculated  in 
the  first  iteration  of  BFIT. 

NCI  -  number  of  records  in  above  arrays 

EEM1  -  intermediate  value  used  in  calculating  the  slope 

significance  statistic. 

K  -  number  of  the  first  record  used  in  the  fit  of  the 

accerelating  creep  data 

00  -  counter  to  terminate  CFIT 

FF  -  slope  significance  statistic 

DW  -  Durbin  Watson  statistic 

MMM  -  LVDT  counter 

ANSM3  -  Interactive  response  that  indicates  whether  or  not 
to  use  CFIT 

NRF  -  number  of  the  last  record  used  in  the  fit  of  the 
decelerating  creep  data 

INTEG 

Integrates  the  power  law  expressions  for  decelerating 
strain  rate  to  obtain  decelerating  strain  for  each  record. 
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SOLUTION 

INTEG  takes  the  fit  parameters  for  decelerating  strain 
rate  from  the  first  call  to  BFIT,  and  performs  an 
integration  to  obtain  the  decelerating  strain  that 
corresponds  to  the  best  fit  strain  rate,  according  to 
equations  [2]  and  [4]  on  page  520  of  Cruden,  1971b. 

The  strain  due  to  accelerating  creep  is  then  obtained 
by  subtracting  the  decelerating  strains  from  the 
experimental  strains. 

This  is  done  for  the  data  from  both  strain  gauges. 

Calling  Sequence:  CALL  INTEG  (El,  E2,  BBO,  BB 1 ,  NRR ,  TT,  L, 
XANS1,  IANS 3 ) 

CF  IT 

The  fit  of  a  power  law  to  the  data  that  is  achieved  in 
subroutine  BFIT  is  a  satisfactory  fit,  not  a  best  fit, 
because  the  fitting  process  is  stopped  once  the  minimum 
statistical  requirements  are  satisfied.  The  fits  that  this 
gives  are  usually  close  to  optimum  for  decelerating  creep, 
but  can  be  significantly  different  from  optimum  for 
accelerating  creep. 

In  CFIT  a  number  of  trial  fits  of  the  accelerating 
creep  that  satisfy  the  minimum  statistical  requirements  are 
examined,  to  identify  the  best  fit  of  the  accelerating  creep 
and  the  best  overall  fit. 
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SOLUTION 

CFIT  calculates  the  strain  rates  for  the  overall  fit, 
by  combining  the  fits  of  the  decelerating  and  accelerating 
components.  It  then  calculates  the  ratio,  R,  of  the  square 
of  the  scatter  of  the  data  points  about  their  mean  to  the 
squares  of  the  scatter  about  this  overall  fit,  a  measure  of 
the  goodness  of  fit. 

The  range  of  data  used  to  calculate  the  accelerating 
creep  fit  is  then  decreased  to  the  next  smallest  range  that 
meets  the  statistical  requirements.  Control  is  then 
transferred  back  to  subroutine  BFIT,  which  fits  a  power  law 
to  this  new  range  of  data.  The  resulting  accelerating  fit 
parameters  are  transferred  to  CFIT,  which  evaluates  the 
overall  fit  obtained  with  the  new  parameters. 

This  process  is  repeated  until  the  overlap  between  the 
accelerating  and  decelerating  creep  goes  to  zero,  or  until 
the  accelerating  strain  rate  come  to  within  one  third  of  the 
decelerating  strain  rate  at  the  beginning  of  the  range  of 
data  used  for  the  fit  of  accelerating  creep. 

The  parameter  and  statistics  corresponding  to  each 
trial  fit  are  tabulated,  to  allow  the  user  to  identify  the 
best  fit. 

Calling  Sequence;  CALL  CFIT  (EEM1 ,  LT,  LE,  W1 ,  BBO,  BB 1 , 

NCI,  K,  00,  FF,  DW,  MMM,  NRR ,  NRF ) 

All  of  these  variables  have  been  defined  previously. 
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CPLOT 

Produces  plots  of  Time  vs  Strain  Rate,  and  Log  Time  vs 
Log  Strain  Rate  with  best  fit  lines. 

SOLUTION 

CPLOT  calculates  strain  rates  for  both  the  experimental 
strain  output  from  CRED,  and  the  accelerating  strain 
calculated  in  INTEG,  and  then  eliminates  all  of  the  negative 
strain  rates. 

If  the  user  desires  a  Time  vs  Strain  Rate  graph,  the 
program  calculates  and  writes  the  extreme  values  for  this 
graph,  and  then  prompts  the  user  for  axis  parameters. 

The  CIVE  Subroutine  GRAPH  is  then  called  to  plot  Time 
vs  Experimental  Strain  Rate. 

If  the  user  desires  a  Log  Time  vs  Log  Strain  Rate 
graph,  the  coordinates  of  the  best  fit  lines  are  calculated, 
and  GRAPH  is  called  to  plot:  experimental  strain  rates, 
deceleration  best  fit  line  and,  if  available,  acceleration 
best  fit  line. 

The  CIVE  subroutine  GRAPH  has  an  interactive  prompt 
sequence  consisting  of  the  following  six  Menu  Options: 

1 .  Plot 

2.  Blow-up 

3.  Redraw 

4.  Subpictures 

5.  Mts-sds 


6. 


Continue 
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The  user  simply  responds  to  the  prompt  by  inputing  the 
number  of  the  desired  option.  Normally  the  user  would  choose 
Option  1,  followed  by  Option  6,  if  a  hard  copy  of  the  plot 
is  desired.  For  more  detailed  information,  consult  the 
documentation  for  GRAPH,  in  the  Library  file: 

CIVE:  GRAPH . DOC . 

One  set  of  graphs  is  produced  for  each  strain  gauge. 

Calling  Sequence:  CALL  CPLOT  (ESTRN1,  ESTRN2 ,  ELPST ,  BBO, 

BB 1 ,  NRR ,  I  SAMI ,  ISAM2,  XANS 1 ,  IANS3,  ANSM1 ) 

Variables  not  defined  previously: 

ANSM1  -  user  response  to  prompts  indicates  if  user  wants 

decelerating  creep  data  plotted  separately. 


The  command  to  execute  this  program  is  as  follows: 

$run  CIVE : CPACK  4=DATA  7=-FILE  T=4SEC 
Note :  A  time  limit  should  be  specified,  and  four  seconds  is 
adequate  for  data  files  with  less  than  500  records. 

The  command  to  obtain  a  hard  copy  of  the  graphs 
generated  in  the  program  is  as  follows: 

$run  *CALCOMPQ  PAR=FI LE=-PDF  FORM=WDWH 
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CPACK  Program  Listing 


1 

2 

3 

4 

5 

e 

7 

8 

9 

10 

1  1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


ccccc 

ccccc 

ccccc 

ccccc 

ccccc 

ccccc 

ccccc 

ccccc 


PROGRAM  CPACK 

A  PACKAGE  OF  PROGRAMS  THAT  WILL  REDUCE 
CREEP  DATA.  AND  FIT  A  POWER  LAW  TO  IT, 
SEPARATING  THE  DATA  INTO  ACCELERATING 
AND  DECELERATING  CREEP 


CCCCC 

CCCCC 

CCCCC 

CCCCC 

CCCCC 

CCCCC 

ccccc 

ccccc 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


DIMENSION  N( 1000) . BBO ( 2 . 2 ) , BB 1 ( 2 . 2 ) . E 1 ( 1000) .E2( 1000) . 

*  TT ( lOOO) , ELPST ( 1000) ,ESTRN1( 1000) ,ESTRN2( 1000) , 

*  LT ( 2 , 1000) , LE ( 2 , 1000) .W1(2, lOOO) . EEM 1 ( 2 ) , NC 1 ( 2 ) 
C 

CALL  CRED  ( DI AM , XLEN . NRR . E 1 . E2 , TT , XANS 1 . I ANS3 , 

*  ITN1 , ITN2, ELPST, ESTRN1 .ESTRN2) 

C 

REAL  Y/'Y'/ 

DATA  MMM/1/ 

00  =  0 
K  =  0 
L  =  1 

NF  I T  =  0 
WRITE ( 7 , 100) 

100  F0RMAT(//20( '  * ' ) . '  FIT  OF  DECELERATING  CREEP  DATA  TO 
•'POWER  LAW  ' , 20( '• ' ) ) 

C 

C  INTERACTIVE  PROMPT 
C 

WR I TE ( 6 , 200 ) 

200  FORMAT* 'WOULD  YOU  LIKE  THE  DECELERATING  CREEP  DATA'/ 
•'PLOTTED  BEFORE  THE  PROGRAM  ATTEMPTS  TO  ISOLATE'/ 
•'ACCELERATING  CREEP?') 

RE AD ( 5 , 300 )  ANSM1 
300  FORMAT ( A  1  ) 

ANSM2  =0 . 

C 

WRITE ( 6 , 202 ) 

202  FORMAT* 'WOULD  YOU  LIKE  THE  PROGRAM  TO  EXAMINE  THE  '/ 
•'ACCELERATING  DATA  FOR  AN  OPTIMUM  FIT?') 

RE AD ( 5 , 302 )  ANSM3 
302  FORMAT* A1) 

C 

CALL  BFIT  (NRR , E 1 . E2 ,TT, L ,N, BBO. BB 1 .XANS1 , IANS3.LT ,LE . 

*  W 1 , E  EM  1 , NC 1 ,K,00,FF,DW, MMM , ANSM3 , NRF , NF I T ) 

C 

IF*  ANSM 1  . EO . Y )  GO  TO  60 
C 

50  CALL  INTEG(E1 ,E2,BB0,BB1 , NRR .TT.L.XANS1 , I ANS3 ) 

C 

L  =2 

WR I TE ( 7 , 101 ) 

101  FORMAT  *//20('*'),'  FIT  OF  ACCELERATING  CREEP  DATA  TO  ', 

• 'POWER  LAW  ' ,20* ' • ' ) ) 

C 

55  CALL  BFIT  ( NRR . E 1 . E2 . TT . L . N , BBO . BB 1 . XANS 1 . I ANS3 . LT , L E , 

*  W1 . EEM1 ,NC 1 ,K,00, FF .DW.MMM, ANSM3.NRF , NFIT ) 


>**  ' 
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. 


61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

B9 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1  10 

1  1  1 

1  12 

1  13 

1  14 

1  15 

1  16 

117 

1  18 

119 

120 
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C 

I F ( ANSM3 . NE . Y )  GO  TO  60 
I F ( 00 . EO . 2 )  GO  TO  60 

CALL  CF IT( EEM 1 , LT , L E , W 1 . BBO , BB 1 . NC 1 , K . 00 . FF , DW , MMM , NRR . NRF ) 
C 

K  =  K+1 
GO  TO  55 
C 

60  CALL  CPL0T(ESTRN1 ,ESTRN2,ELPST,BB0,BB1 ,  NRR ,  I  TNI  ,  ITN2, 

*  XANS1 . IANS3, ANSM1 . ANSM2) 

C 

I F ( ANSM 1 . NE . Y )  GO  TO  999 
I F ( NF I T  . EO. 1 )  GO  TO  555 
C 

C  INTERACTIVE  PROMPT 
C 

WR I TE ( 6 , 20 1 ) 

201  FORMAT( 'WOULD  YOU  LIKE  THE  PROGRAM  TO  ATTEMPT  TO  FIT'/ 
•'YOUR  DATA  TO  ACCELERATING  CREEP  AS  WELL?') 

READ(5 , 301 )  ANSM2 
301  FORMAT ( A  1 ) 

IF(ANSM2 .NE . Y)  GO  TO  999 
ANSM 1 =0 . 

GO  TO  50 
555  WR I TE ( 6 , 40 1 ) 

401  FORMAT ( ' THE  ENTIRE  SET  OF  DATA  IS  FITTED  TO  THE'/ 
•'DECELERATING  CREEP') 

999  STOP 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


ccccc  ccccc 

CCCCC  SUBROUTINE  CRED  CCCCC 
CCCCC  REDUCES  EXPERIMENTAL  DATA  CCCCC 
CCCCC  CCCCC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  CRED  ( DI AM . XL EN . NRR . E 1 . E2 . TT . XANS 1 . I ANS3 , 

•I  TNI . ITN2, ELPST.ESTRN1 , ESTRN2 ) 

DIMENSION  RLOAD( 1000) , RDEF 1 ( 1000) , RDEF2( 1000) , IMN( 1000) , 

•  TIME ( 1000) , VDEF 1 ( 1000) , VDEF2( 1000) . IDY( 1000) , 

•  E  STRN2 ( 10O0) , RCELLP ( 1000) , ESTRN1 ( 1000) , 

•  IHR ( 1000) , IMI ( 1000) , ISE( 1000) , ELPST( 1000) , 

•  E 1 ( lOOO) , E2( 1000) , TT ( 1000) ,VAR( 12, 1000) , XNAME ( 12) 

C 

C 

REAL  Y/'Y'/.H/'H'/.M/'M'/.TH/'TH.  '/.R1/'R1  ,  '/,R2/'R2,  '/, 

•DU/ 'DU, '/.TM/'TM, '/,TN/'TN, '/.TD/'TD, '/,TS/'TS, '/.FACT/1 . 0/ . 
• Y  Y  Y / ' y ' / , NNN / ' n ' / 

REAL  IMN/ 1000*0 . / , IDY/ 1000*0 . / , IHR/ 1000*0 . / , 

•  I M I / 1 000 *0./, ISE/ 1 0OO  *  0 . / 

C 

C  PROMPTS  FOR  INPUT  PARAMETERS 
C 

WR I TE ( 6 , 150) 

150  FORMAT ( ' THI S  PROGRAM  ACCEPTS  DATA  IN  THE  FOLLOWING  FORMAT:'/ 
•'LABEL, TIME  CLOCK, 3  COUNTERS , LOAD . CELL  PRESSURE ,DEF . 1 ,DEF . 2 ' 
•/'IS  YOUR  DATA  IN  THIS  FORM?  ( Y , N ) ' ) 

11  RE AD ( 5 , 250 ) X ANS 1 
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121 

122 
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124 
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127 

128 

129 

130 

131 
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136 

137 

138 

139 

140 

14  1 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 


non  o  o  o 
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250  FORMAT ( A  1 ) 

IF(XANS1 .NE . YYY. AND.XANS1 . NE .  NNN )  GO  TO  12 
WR I TE ( 6 , 26  1  ) 

261  FORMAT( ' ‘ERROR*  ALL  RESPONSES  MUST  BE  IN  UPPER  CASE  LETTERS.' 
‘/'RE-ENTER  RESPONSE  TO  PREVIOUS  QUESTION:') 

GO  TO  11 
12  WR I TE ( 6 , 151 ) 

151  FORMAT ( ' DO  YOU  HAVE  AN  LVDT  CALIBRATION  FACTOR?') 

RE AD (5, 251 )XANS2 

251  FORMAT ( A  1  ) 

IF(XANS2 .NE . Y)  GO  TO  199 
WR I TE ( 6 , 160) 

160  FORMAT( ' INPUT  THE  LVDT  CALIBRATION  FACTOR  (MULTIPLIER):') 

RE AD( 5 , 260) FACT 
260  FORMAT ( F 1 0 . 5 ) 

199  WR I TE ( 6 , 152) 

152  FORMAT( ' INPUT  DI AMETER ( MM ) , L ENGTH( MM ) . AND  WEIGHT (GM) , ' 
‘/'SEPARATED  BY  COMMAS:') 

READ( 5 , 252 )DIAM , XLEN , WEIG 

252  FORMAT ( 3F 10 . 3 ) 

WR  I TE ( 6 , 153) 

153  FORMAT( ' INPUT  SAMPLE  NUMBER  (8  CHARACTER  MAX.):') 

RE AD (5, 253) I  SAMI . ISAM2 

253  FORMAT ( 2A4  ) 

WRITE ( 6 , 154) 

154  FORMAT( ' INPUT  TEST  NUMBER  (8  CHARACTER  MAX.):') 
READ(5,254)ITN1 .ITN2 

254  FORMAT ( 2A4 ) 

WRITE  INPUT  PARAMETERS 


WR ITE(7,20)ITN1 , ITN2 

20  FORMAT(/, 'TEST  NUMBER 
WR I TE ( 7 , 2 1 )  ISAM1.ISAM2 

21  FORMAT(/ ,' SAMPLE  NUMBER 

WR I TE ( 7 , 22  )  XLEN, DIAM, WEIG 

22  FORMAT(/, 'SAMPLE  LENGTH 

»  /. 'SAMPLE  DIAMETER  = 

»  /. 'SAMPLE  WEIGHT 

I F ( XANS 1 . EQ . Y )  GO  TO  23 


2A4  ) 

2A4  ) 

F10.3, ' 

mm '  , 

F10.3. ' 

mm '  % 

F10.3. ' 

gm '  ) 

PROMPT  FOR  INPUT  FORMAT  (NON-STANDARD) 


WR I TE ( 6 , 155) 

155  FORMAT( ' INPUT  NUMBER  OF  ELEMENTS  IN  EACH  LINE  (+  COMMA):') 

RE  AD ( 5 , 255 ) NV 

255  FORMAT (12) 

WR  I TE ( 6 , 156) 

156  FORMAT( ' INPUT  NUMBER  OF  STRAIN  GAUGES  (LVDTS);(1  OR  2):') 
READ(5,256)IANS3 

256  FORMAT (II) 

WR I TE ( 6 , 157) 

157  F0RMAT(5X , 'THE  ONLY  DATA  THAT  THIS  PROGRAM  REQUIRES  FOR  INPUT' 
‘/'IN  DEVICE  4,  IS  TIME  AND  DISPLACEMENT.  LIST  THE  ELEMENTS  IN'/ 
•'ONE  LINE  OF  YOUR  INPUT  DATA,  IN  PROPER  ORDER,  USING  "R1"'/ 
•'FOR  THE  FIRST  STRAIN  GAUGE  READING.  "R2”  FOR  THE  SECOND  (IF'/ 
♦'THERE  IS  ONE),  AND  "TN" , "TD" , "TH" , "TM" , "TSM ,  FOR  THE  TIME'/ 
•'IN:  MONTHS, DAYS. HOURS. MINUTES. AND  SECONDS  (YOU  MAY  HAVE  ONLY' 
•/'ONE  OF  THESE  TIME  PARAMETERS).  USE  THE  DUMMY  VARIABLE  "DU"'/ 
•'FOR  ALL  OTHER  VALUES.  (SEPARATE  THE  ELEMENTS  WITH  COMMAS)') 
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181 
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192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

210 
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212 

213 

214 

215 

216 

217 

218 

2  19 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 


non  no  non 
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READS  NON-STANDARD  DATA 


257 

5 


8 

9 

16 

6 

17 
7 

18 

31 

42 

32 

43 

33 

44 

34 

45 

35 

46 
39 

300 


READ(5,257)(XNAME(J) . J=1 ,NV) 

FORMAT ( 12A3) 

L  =0 
L  =  L+  1 

RE AD ( 4 , 300 . END =8 )(VAR(K,L) ,K=1 .NV) 
N  =  L 


GO  TO  5 
K  =  0 


K  =  K+1 


I F ( XNAME (K ) .NE .TH) 
DO  16  L* 1 , N 
I  HR ( L ) * VAR ( K , L  ) 

GO  TO  39 

I F ( XNAME ( K ) .NE .R1 ) 
DO  17  L  = 1 , N 
RDEF1(L)=VAR(K.L) 

GO  TO  39 

I  F ( XNAME ( K ) .NE ,R2) 
DO  18  L  =  1  ,  N 
RDEF2(L)=VAR(K,L) 

GO  TO  39 

I F ( XNAME ( K ) .NE .TN) 
DO  42  L  = 1 , N 
I MN( L)=VAR(K,L) 

GO  TO  39 

I F ( XNAME ( K ) .NE . TD ) 
DO  43  L  =  1 . N 
IDY ( L ) =  VAR ( K , L ) 

GO  TO  39 

I F ( XNAME ( K ) .NE .TH) 
DO  44  1= 1 , N 
I HR (L)=VAR(K,L) 

GO  TO  39 

I F ( XNAME ( K ) .NE . TM ) 
DO  45  L  = 1 , N 
IMI(L)=VAR(K,L) 

GO  TO  39 

I F ( XNAME ( K ) .NE .TS) 
DO  46  L  = 1 , N 
I  SE ( L ) = VAR ( K , L  ) 

IF(K . LT .NV)  GO  TO  9 
FORMAT ( 12G20) 


GO  TO  6 


GO  TO  7 


GO  TO  31 


GO  TO  32 


GO  TO  33 


GO  TO  34 


GO  TO  35 


GO  TO  39 


56  DO  51  K= 1 , N 
RLOAD ( K ) =0 . 0 

51  RCELLP ( K ) =0 . 0 

I F ( I ANS3 . EO . 2 )  GO  TO  55 
DO  52  K= 1 , N 

52  RDE F 2 ( K ) =0 . 0 
GO  TO  55 

READS  STANDARD  DATA 


23  J  =  0 

WRITE ( 6 , 158) 

158  FORMAT ( 'DOES  YOUR  TIME  CLOCK  HAVE  A  DOUBLE  DIGIT  MONTH?') 
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262 

263 

264 

265 

266 

267 

268 
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270 

271 

272 
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274 

275 

276 

277 

278 

279 

280 

281 

282 
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READ(5.258)XANS3 
258  FORMAT ( A  1 ) 

30  d  =  d+1 

I F ( XANS3 . EQ . Y ) GO  TO  47 

READ(4,40, END =55 )  IMN( d ) . IDY ( d ) , IHR( d ) , IMI ( d ) . I SE ( d) . ICTR. 

*  II, IO.RLOAD(d) .RCELLP(d) , RDEF 1(d) ,RDEF2(d) 

40  FORMAT (7X.5G2. IX. I5.2I3.4G20) 

GO  TO  48 

47  READ(4,49, END =55 )  IMN( d) . IDY ( d ) , IHR ( d) , IMI ( d ) . I SE ( d ) , ICTR . 

*  II , IO.RLOAD(d) , RCELLP ( d ) .RDEF 1 (d) ,RDEF2(d) 

49  FORMAT (8X.5G2, IX. 15,213, 4G20) 

48  N  =  d 

GO  TO  30 


55  RILOAD  =  RLOAD(  1  ) 

RCELLP ( 1 )=RCELLP(  1  ) 

ALVDT  1 =RDEF  1  (  1  ) 

BLVDT2=RDEF2(  1  ) 

WRITES  INITIAL  READINGS 

WRITE (7, 70)  RILOAD, RCELLP ( 1 ) , ALVDT 1 .BLVDT2 
70  FORMAT(/, ' INITIAL  LOAD  =  '.F10.3.'  KN ' . 

*  /.'INITIAL  CELL  PRESSURE  =  '.F10.3,'  KPA ' , 

*  /.'INITIAL  READING  FOR  LVDT1  =  '.F17.10.'  mm', 

*  /.'INITIAL  READING  FOR  LVDT2  =  '.F17.10,'  mm') 

CONVERTS  TIME  FROM  CLOCK  READING  TO  HOURS 

59  IMNST  =  IMN( 1 ) 

DO  50  d= 1 ,N 

60  TIME(d)=(IDY(d)*24. )+ IHR ( d ) ♦ ( ( IMI ( d ) +1 SE ( d ) /60 . )/60. ) 

I F ( IMN(d) . EO . IMN( 1 ) )  GO  TO  50 

I F ( IMN( 1 ) .NE . 4 . OR . IMN( 1 ) . NE .6 . OR . IMN( 1 ) .NE .9 .OR . IMN( 1 ) .NE . 1 1 ) 
•GO  TO  61 
IDY(d)=IDY(d)+30 
GO  TO  63 

61  I F ( I MN( 1 ) . NE . 2 ) GO  TO  62 
IDY(d)=IDY(d)+28 

GO  TO  63 

62  IDY(d)=IDY(d)+31 

I F ( I MN ( 1) .EO. 1 2 ) GO  TO  64 

63  I MN ( d )  =  I MN (  d ) -  1 
GO  TO  60 

64  I MN ( d ) = 1 2 
GO  TO  60 

50  CONTINUE 

INITIALIZES  VARIABLES 

ELPST ( 1 )=0. 

VDE  F 1 ( 1)=0. 

VD  E  F  2 ( 1)=0. 

ESTRN1 ( 1)=0. 

ESTRN2 ( 1)=0. 

CALCULATES  ELAPSED  TIME.  LOAD,  CELL  P..  AND  ENGG.  STRAIN 
DO  80  d=2 , N 
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342 
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VDEF 1 ( J) =RDEF1 ( J)-RDEF 1(1) 

VDE  F  2(J)=RDEF2(J)-RDEF2( 1 ) 

ELPST( J)=TIME( J)-TIME( 1 ) 

RL0AD( J)=RL0AD( J)-RIL0AD 
RCELLP(U)=RCELLP( J)-RCELLP(  1  ) 

ESTRN1 ( J)=( VDEF 1 (U)/XLEN) *1000000. * FACT 
ESTRN2 (U)=(VDEF2(J)/XLEN)* lOOOOOO . ‘FACT 
80  CONTINUE 
RLOAD ( 1  )=0. 

C 

C  WRITES  TIME,  LOAD.  CELL  P.,  AND  ENGG.  STRAIN 
C 

WR I TE ( 7 , 90 ) 

90  FORMAT(/, IX. 'NUMBER' ,4X, 'TIME ' ,9X, 'LOAD' ,8X, 'CELL  PRESSURE'. 

*  2 ( 3X . ' ENGG .  STRAIN  ')./. 1 IX ,' HOURS '. 9X ,' KN '. 14X , 

*  ' KP A ' , 9X , ' # 1  (MICRO) ' ,7X, ' H2  (MICRO)') 

DO  100  J=1 , N 

WRITE (7, 1 10)  U . ELPST ( J ) , RLOAD( J ) , RCELLP ( J ) , ESTRN1 ( J) , ESTRN2 ( J ) 
110  F0RMAT(2X,I3.3X,F9.3,3X.E13.6.2X.E13.6,3X.F11 . 3.6X.F11 .3) 
lOO  CONTINUE 
C 

C  DELETES  FIRST  LINE  OF  DATA 
C 

NRR  =  0 

NRR=N- 1 

DO  85  1=1 , NRR 

E  1  (  I  )  =ESTRN1 ( 1+ 1 ) 

E2(I ) =E  STRN2 (1  +  1 ) 

85  TT ( I ) =ELPST (1+1 ) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


CCCCC  CCCCC 
CCCCC  SUBROUTINE  BFIT  CCCCC 
CCCCC  FITS  A  POWER  LAW  TO  THE  DATA  CCCCC 
CCCCC  BY  RESTRICTION  ON  DW  S  FF  CCCCC 
CCCCC  CCCCC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE  BFIT  ( NRR . E 1 . E2 . TT , L . N . BBO , BB 1 , XANS 1 , I ANS3 , LT , L E . 

*  W1 , EEM1 .NCI ,K,00, FF .DW.MMM. ANSM3.NRF ,NFIT ) 
DIMENSION  T ( 1000) . A( 10) , AT( 1000) , E 1 ( 1000) . E2( 1000) . 

*W( 1000) . B ( lOOO) , E A (  lOOO) , AE ( 1000) ,N( 1 ) , BBO (2. 2) ,BB1 (2.2)  . 
*DENT ( 48 ) , VDE F ( 1000) , DLOAD ( lOOO) ,TT( lOOO) , DUMMY ( 1000)  . 

*  E (  1000) . E  E ( 1000) . ER ( 1000) .LT(2,  1000) ,LE(2. lOOO) ,W1 (2. 1000) , 
*EEM 1(2) , NC 1(2) , NRF ( 2 ) 

REAL  Y/'Y'/, BLANK/'  '/.LT.LE 
C 
C 

Z=0 

IF(OO.GT.O)  Z=1 
IF(OO.EQ.O)  K= 1 
MM  =  0 

700  MM=MM+ 1 

IF(ANSM3.E0. Y.AND.L.E0.2)  MM=MMM 
NR=NRR 

I F ( MM . EO . 1 ) GO  TO  710 
DO  705  11=1, NR 


■ 


*“T  3T  ■  A  t'l  * 


>  03 


.  -  *i 


- 


361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

37  1 

372 

373 

374 

375 

376 

377 

378 

379 

380 

38  1 

382 

383 
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387 

388 

389 

390 

391 
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393 
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non  ono  ooo 


113 


is»(  1 1 )  =  1 1 

705  DUMMY ( II )=E2( II ) 

GO  TO  720 
710  DO  715  11=1 .NR 
M( II  )  =  II 

715  DUMMY (II)=E1(II) 

720  CONTINUE 

PROMPT  TO  DETERMINE  LIMIT  FDR  SLOPE  SIGNIFICANCE 

I F ( L . E  0 ■ 2 . OR . MM . E  0 . 2 )  GO  TO  515 
WRITE (6 , 200) 

200  FORMAT( ' INPUT  DESIRED  LIMIT  FOR  TEST  OF  SLOPE  SIGNIFICANCE  IF' 
♦/'OTHER  THAN  10.  (REAL  NUMBER,  TERMINATED  WITH  A  COMMA):') 
READ( 5 , 100) FFL 
100  FORMAT ( G4 ) 

I  F ( FFL .LT .0. 1 )  FFL= 10 . 

SMOOTHS  DATA  TO  INCREASING  POSITIVELY 

515  I =K- 1 
J  =  K-  1 
12  1=1+1 
d  =  J+1 

AE ( I  )=DUMMY(J) 

AT(I )=TT(J)*60. 

W( I  )  =  1 .0 

11  I F ( I  . EQ  K )  GO  TO  12 
14  IF(AE(I )— AE( 1—1 ) )  13,  13.  4 

13  AE( 1-1 ) = ( AE ( I )*W( I )+AE(I-1)*W( 1-1 ) )/(W( I )+W(I-1 ) ) 

AT(I-1)=(AT( I )*W( I )+AT( 1-1 )*W( 1-1) )/(W( I )+W(I-1 ) ) 

W(I-I) =W( I )+W( 1-1 ) 

1  =  1-1 

IF(I-K-I)  12,14,14 
4  CONTINUE 

I F ( J . LT . NR ) GO  TO  12 

INITIALIZES  VARIABLES 


ww=o . 

BB  =  0 . 
CONB 1 =0 . 
C0NB0=0 . 
TE  =0 . 

DW  =  0 . 
EER=0 . 
EEM=0. 
EES=0 . 
SUMT=0 . 
SUMET  =0 . 
SUMT2=0. 
SUME  =0 . 
SUME  2  =  0 . 
SXX=0.0 
DWW=0 . 

SX  Y  =0 . 0 
WWA  =  0 . 

AF  =0 . 

M  =  K+  1 
NC  =  I 
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CALCULATES  LOG  STRAIN  RATE,  LOG  TIME 
DO  1  J=M,NC 

ER(U)=(AE( J)-AE(U-1))/(AT( J)-AT( J-1 ) ) 
T(J)=(AL0G10(AT(J) )+AL0G10( AT ( J- 1 ) ) )/2.0 

6  W(J)=(W(J)+W(J-1))/2. 

E( J)=AL0G1O(ER( J) ) 

WW=WW+W( J) 

EE ( J) =0 . 

SUMT  =  SUMT  +T(J)*W(J) 

SUME=SUME+E( J)*W(d) 

SUMET  =SUMET+E(J)*T(J)*W(J) 

SUME2  =  SUME2-*-E  (J)*E(J)*W(J) 

1  SUMT 2= SUMT 2 +T  (J)*T(J)*W(J) 
SUME2=SUME2-SUME*SUME/WW 
SUMT  2  =  SUMT  2 - SUMT  *  SUMT / WW 
SUMET  =SUMET-SUME*SUMT/WW 
FME  =  SUME/WW 
FMT=SUMT/WW 
DO  7  J=M,NC 

SXY=SXY+(T( J)-FMT)*(E( J)-FME)*W( J) 

7  SXX=SXX+W( J)*(T( J)-FMT)**2 


USES  LEAST  SQUARES  TO  FIT  DATA 


B 1 =SXY/SXX 

BO=FME-FMT*B 

1 

DO  2  J=M , NC 

EE  ( <J  )  =BO+B  1  * 

T(  J) 

EES=EES+( EE( 

U)-E( J) 

) *W( J ) 

EA( J)=E( J)-E 

E(U) 

EER=EER+W( J) 

*  (  E  E  (  J ) 

-E( J) )**2 

EEM=EEM+W( J) 

* ( EE ( J) 

-FME ) *  *2 

I F  (  J . L  E . M ) GO 

TO  2 

DW=DW+( EA ( J) 

-EA( J-1 

)  )**2 

DWW=DWW+E A ( J 

)  *  E  A  (  J  ) 

2  CONTINUE 
SSDYX  =  E  E  R 
EER=EER/(WW-2 . ) 

FF=EEM/EER 

CONB 1 =SQRT( EER/SUMT2 ) 

CONBO=CONB 1 *SQRT ( ( SUMT2* WW+5UMT*SUMT ) / ( WW* WW ) ) 
DW=DW/DWW 

TRANSFER  VARIABLES  FOR  SUBROUTINE  CFIT 

IF(00.EQ.2. AND . Z . EQ . 0)  00=0 

IF (NR  NE .NRR . OR . L .NE . 1 )  GO  TO  111 

E  EM  1 ( MM ) =EEM 

NC 1 ( MM ) =NC 

DO  21  J=1 ,NC 

LE(MM,J)=E(J) 

LT(MM,U)=T( J) 

21  W1(MM, J)=W( J) 

CHECKS  DW  AND  FF 


111  NDF  =  NR -K+ 1 
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IF(NDF . LE . 20)  DU  = 1 . 36+ ( NDF -  1 5 ) *0 . 01 

I F ( NDF . GT . 20. AND .NDF .LE . 30)  DU= 1 . 4 1 + ( NDF -20 ) *0 . 008 
IF(NDF .GT .30. AND.NDF.LE .40)  DU=1 . 49+ ( NDF -30 ) *0 . 005 
IF (NDF .GT .40. AND  NDF . LE . 50)  DU= 1 . 54+ ( NDF -40) *0.005 
IF (NDF .GT . 50. AND .NDF . LE .60)  DU= 1 .59+ (NDF -50) *0.003 
IF (NDF .GT .60. AND .NDF .LE .95)  DU=  1 . 62+ ( NDF -60) *0 . 002 
IF (NDF  . GT.95)  DU= 1 .69 
IF(DW.GT.DU.AND.FF.GT.FFL)  GO  TO  817 
I F ( L . EO -  1)  NR=NR- 1 
I F ( L . EO . 2 )  K=K+ 1 
GO  TO  515 


NDI FF  =NRR-NR 
IF(NDIFF . EO.O)  NFIT=1 
I F ( ANSM3 . NE . Y)  GO  TO  8 
IF(L . EO. 1 )  GO  TO  8 
I F ( 00 . EO -  2 )  GO  TO  8 
BBO ( L , MM ) =B0 
BB 1 ( L , MM ) =B 1 
GO  TO  702 

WRITES  DATA  AND  STATISTICS 

8  WR I T  E ( 7 , 206 ) MM 

206  F ORMAT ( ///30X , '  DATA  FROM  LVDT  NO. ',12,//) 

416  CONTINUE 

WR I TE ( 7 , 208 ) 

208  F0RMAT(9X , 'TRANSFORMED  DATA') 

WRITE ( 7 , 209 ) 

209  FORMAT( 1 IX, 'TIME ' ,5X, 'STR  RATE.  E'.5X.'L0G  E',9X,'L0G  EE'.4X, 
* 'LOG  E  -  LOG  EE ' ,4X, 'W' ) 

WR ITE(7,216) 

216  FORMAT ( 1 1 X , ' ( MI N ) ' , 3X , ' (MICRO . E/MIN) ' ) 

DO  9  J=M , NC 
T ( J ) *  10*  * ( T ( J ) ) 

WR I TE ( 7  ,  106 )T(U) ,ER(J) ,E(J) ,EE( J) , E A ( J )  ,W(J) 

106  F0RMAT(7X,F 12 .4, 2X,F 12 .6,3( 1X,E13.6),F7.2) 

9  CONTINUE 

WR ITE(7,210) 

210  FORMAT( // , 9X , ' FIT  PARAMETERS') 

WR I TE ( 7 , 105 )B0 

105  FORMAT ( 1 1 X ,  '  INTERCEPT,  BO '  ,  1 6X . E  1  5 . 6 ) 

WR I TE ( 7 ,  1  15)B  1 
WR I TE ( 7 , 1 58 ) CONBO 
WRI TE ( 7 , 157) CDNB 1 
WRITE  ( 7 , 159 ) DW 


157 

FORMAT ( 1 IX , ' 

CONFIDENCE  LIMIT  ON  B1  T. 

'  ,F12.7) 

115 

FORMAT (  1  IX  ,  ' 

SLOPE.  B 1 

'  .F12.7) 

158 

FORMAT (  1  IX,  ' 

CONFIDENCE  LIMIT  ON  BO  T. 

'  ,F  12 .7) 

159 

FORMAT (  1  IX .  ' 

DURBIN  WATSON  STATISTIC 

'  , F7 . 3 ) 

WR I T  E ( 7 ,  180) 

NDF 

1  BO 

FORMAT ( 1  IX ,  ' 

DEGREE  OF  FREEDOM  FOR  DW 

'  ,  13) 

WR I TE ( 7 ,  181  ) 

DU 

181 

FORMAT (  1  IX,  ' 

UPPER  LIMIT  FOR  DW 

'  ,  F  7 . 3  ) 

WR I TE ( 7 , 156) 

FF 

156  FORMAT ( 1 1 X , '  TEST  OF  SLOPE  SIGNIFICANCE  '.F9.3) 
WRITE  (7.211) 

211  F0RMAT(//,9X, 'DATA  FOR  COMPARISON  TESTS') 
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WW= WW-2 . 

WR I TE ( 7 , 122) WW 

122  FORMAT ( 1 1 X , '  WEIGHTING  '.F6.1) 

WR I T  E ( 7 , 104 ) FME , FMT 
WRITE (7 , 214 )SUME2 , SUMT2 

104  FORMAT ( 1 1 X , '  MEAN  STRAIN  '.F11.3. 

*'  MEAN  TIME  '  ,F11 .3) 

WR ITE(7,213) SUMET , SSD YX 

214  FORMAT ( 1 1 X , '  SSDY  '.5X.F11.3,'  SSDX  '.F11.3) 

213  FORMAT ( 1 1 X , '  SPDXY ' . 5X , F 1 1 . 3 . '  SSDYX  '.F11.3) 

WRITE(7 ,212) 

212  F0RMAT(/,9X. 'CHECK' ) 

WR I TE ( 7 , 107) EES 

107  FORMAT ( 1 1 X , '  SUM  OF  RE S IDUALS ' . 14X . F 1 2 . 6 ) 

WRITE (7 , 124  )  K 

124  FORMAT ( /// , 1  1  X  ,  '  DATA  STARTS  AT  #',I3) 

WRITE ( 7 , 125)  NR 

125  FORMAT ( 1 IX , '  DATA  ENDS  AT  #',I3) 

C 

C  INITIALIZES  COUNTERS  FOR  SUBROUTINE  CFIT,  IF  USED 
C 

I F ( L . EO . 1 )NRF ( MM ) =NR 
IF(ANSM3.EQ.Y. AND . L . EQ . 2 )  GO  TO  704 
GO  TO  703 

704  MMM=2 
K  =  1 
Z=0 

c 

C  DIRECTS  PROGRAM  ACCORDING  TO  PROGRAM  OPTIONS  IN  EFFECT 
C 

703  BBO ( L , MM ) =B0 
BB 1 ( L , MM ) =B 1 
IF (XANS1 . EO. Y)  GO  TO  701 
IF( IANS3 . EO . 1 )  GO  TO  702 

701  I F ( MM . EO . 2 ) GO  TO  702 
GO  TO  700 

702  CONTINUE 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


CCCCC  CCCCC 
CCCCC  SUBROUTINE  INTEG  CCCCC 
CCCCC  INTEGRATES  STRAIN  RATE  TO  DETERMINE  CCCCC 
CCCCC  ACCELERATING  STRAIN  CCCCC 
CCCCC  CCCCC 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  INTEG( E 1 ,E2,BB0,BB1 , NRR , TT , U J , XANS 1 , IANS3) 
DIMENSION  EDI FF ( 1000) ,E( 1000) . EE ( 1000) ,E1( 1000) , BBO (2, 2) , 

*  BB1(2.2) ,TT( lOOO) ,E2( 1000) 

REAL  Y/'Y'/ 

C 

c 

MM  =  0 
NR=NRR 
20  MM=MM+ 1 

I F ( MM . EO . 2 )  GO  TO  43 
DO  40  1=1 , NR 
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40  e ( i )  =  e 1 1 1 ) 

c 

C  INTEGRATES  STRAIN  RATE  EQUATION.  TO  GET  STRAINS 
C 

43  BO=BBO( JJ . MM ) 

B  1  =BB 1 ( JJ . MM ) 

B0= 10.  *  *  BO 
B 1 =B  1+  1 
BO  =  BO/B  1 
I F ( B 1  )  12.  12,  13 

12  DO  9  I = 1 .NR 

T  E  R  M  =  0 . 00000 1 / T  T ( I  ) 

9  EE(I)=-(B0)M10**(-6*B1))*( 1-(TERM**(-B1))) 

GO  TO  2 

13  DO  1  I = 1 . NR 

1  EE(  I  )=B0*( (TT( I )**B1  )-(  1 0*  * ( - 6  *B  1  )  )  ) 

C 

C  SUBTRACTS  DECELERATING  STRAIN  FROM  EXPERIMENTAL  STRAIN 
C 


c 

c 


2  DO  21  I = 1 .NR 
2  1  ED  I F  F ( I  )  =E(  I  )-EE(  I  ) 


30 


15 
3  1 
999 


I F ( MM . E Q . 2 )  GO  TO  15 
DO  30  1=1. NR 
E  1  (  I  )  =  ED I F  F (  I ) 

E( I ) =E2( I  ) 

I F ( X ANS 1  . E Q . Y )  GO  TO 
I F ( IANS3 . EQ. 1 )  GO  TO 
GO  TO  20 
DO  31  I = 1 . NR 
E2 ( I ) =ED I F  F ( I ) 

RETURN 

END 


20 

999 
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SUBROUTINE  CFIT 
OPTIMIZES  THE  FIT  OF  THE 
POWER  LAW  EXPRESSION  TO  THE 
ENTIRE  RANGE  OF  DATA 


CCCCC 

CCCCC 

CCCCC 

CCCCC 

CCCCC 

ccccc 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


SUBROUTINE  CFIT(EEM1.LT.LE.W1  .BB0.BB1  , NC 1 . K . 00 . 

*  F  F , DW . MMM . NRR . NRF ) 

DIMENSION  L T ( 2 .  1000) ,LE(2.  1000)  .W1(2.  1000) . 

*  BBO ( 2 ,2).BB1(2.2).EEM1(2), NC 1(2), NRF ( 2 ) 

REAL  LT.LE 

C 

C 

MM=MMM 
WW  =  0 
M  =  0 
M  =  K+  1 

NC  =NC 1 (MM) 

E  E  R  =0 
C 

C  CALCULATES  R  FOR  CURRENT  ACCELERATING  CREEP  PARAMETERS 
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c 

DO  1  J  =  2  .  NIC 

DE  =BBO  (  1  ,MMMUT(MM.  J)*BB1<  1  .  MM  )  ) 

A  E  =  BBO ( 2 , MM  )  +  ( L  T ( MM. J ) *BB 1 ( 2, MM) ) 

E  E  =  ALOG  1 0  (  (  10*  *  AE  )  ■*•  (  10*  *DE  )  ) 

EER=EER+W1(MM.J)*(EE-LE(MM.J))*‘2 
1  WW=WW* W 1 ( MM , J ) 

EER=EER/(  WW-  2  .  ) 

R=EEM1(MM)/EER 
I F ( 00 . NE .0)  GO  TO  2 
C 

C  WRITES  TABLE  OF  FIT  STATISTICS 
C 

WRI TE ( 7 . 1 OO ) MM 

100  FORMAT (//.22X,  'PARAMETERS  FOR  EVALUATING  TOTAL  FIT.'. 

*  '  LVDT  NO  .  :  '.II.// 

*34X. 'ACCELERATING  CRE E P ' . 1 8X . ' TOT AL  FIT'/ 

*16X.  'RANGE  INTERCEPT, BO  SL OP E . B  1  '  . 6 X .  ' R '  . 8 X .  ' DW '  . 8 X .  ' R  '  ) 

2  WRITE(7, 10 1 ) K , NRR . BBO( 2 . MM ) . BB 1 ( 2 . MM ) . F F . DW . R 

101  FORMAT (  1 4X ,  I  3 ,  -  '  . I  3 . E 15 . 6 , F 12 . 7 , 2F9 . 3 . 2X . F9  3) 

C 

C  TEST  TO  SEE  IF  DATA  MEETS  TERMINATING  CRITERIA 
C 

EA=BBO( 2 ,MM)+( LT(MM.M) ‘BB 1 ( 2 . MM ) ) 

ED=BBO ( 1 .MM)*(LT(MM.M)*BB1( 1 , MM ) ) 

E  A  =  E  A+0 .  1 
00=  1 

I F ( EA . GE . ED . OR . M . GE . NRF (MM) )  GO  TO  3 
GO  TO  4 

3  00  =  2 
K  =  K-  1 

4  RETURN 
END 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


ccccc 

CCCCC 

ccccc 

SUBROUTINE 

CPLOT 

CCCCC 

ccccc 

MAKES  PLOTS 

OF  : 

CCCCC 

ccccc 

STRAIN  RATE 

VS  TIME 

CCCCC 

ccccc 

AND 

CCCCC 

ccccc 

ccccc 

LOG  STRAIN  RATE 

VS  LOG  TIME 

CCCCC 

ccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE  CPLOT ( ESTRN 1 . ESTRN2 . ELPST . BBO . BB 1 . NRR . 

*  I  TNI . ITN2.XANS1 . IANS 3. ANSM1 . ANSM2 ) 

DIMENSION  ESTRN 1(1), ESTRN2 (  1  )  . ELPST (  1  ) . BBO( 2 . 2 ) . BB  1  ( 2 . 2  )  . 

‘  A  E (  1000) . E (  1000) , T (  10OO) ,AT(  1000) , DX ( 4 ) . D Y ( 4 ) . 

*  LABELS(24) ,OPTNS( 25) ,TE( 10O0) 

REAL  LOG/'LOG  ' / , Y/ ' Y ' / . TWO/ ' TWO  ' / . ONE / ' ONE  '/.BLANK/'  '/ 

c 

C  INITIALIZE  VECTORS 

C 

REAL  LABELS/'TEST' . '  NO.'.':  '.2*'  '.'  LVD'.'T  NO'. 

*'•:  '.4*'  '.'TIME'. 3*'  '.'STRA'.'IN  R'. 

• 'ATE  ' .4* '  '/. I TN 1 . I TN2 

REAL  OPTNS  / 1  . O . 24  * ' NO  '/ 

c 

C  DETERMINE  IF  IN  BATCH  MODE 
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CALL  CREPLY(&8) 

GO  TO  9 

8  OPTNS ( 4 ) = 1 
C 

C  PROMPTS  TO  DETERMINE  DESIRED  PLOTS 

C 

9  WR I T E ( 6 , 200 ) 

200  F ORMAT (  WOULD  YOU  LIKE  A  LOG  STRAIN  RATE  -  LOG  TIME  PLOT?(Y.N)') 

READ( 5 . 300) ANS 1 

300  FORMAT ( A  1  ) 

WR I T  E ( 6 , 20 1 ) 

201  FORMAT ( '  WOULD  YOU  LIKE  A  STRAIN  RATE  -  TIME  PLOT?(Y.N)') 

READ( 5 , 301 ) ANS2 

301  F  ORMAT ( A  1 ) 

I F ( ANS  1  . NE  . Y . AND . ANS2 . NE  . Y )  GO  TO  1000 
I F ( ANSM2 . EO . Y )  GO  TO  111 
C 

C  LOAD  GRAPH 
C 

INTEGER  LOADF ,  GRAPH,  L SWS/Z00800040/ 

GRAPH  =  LOADF  (  ' C I VE : GR A PH+ *  I G+ * PLOT L I B  O,  LSWS,  O  ) 

c 

c 

c 

1 1 1  MM  =  0 
50  MM  =  MM*  1 

I F ( MM . EO . 2 ) GO  TO  1 
LABE  L  S ( 4 )  =  I TN 1 
LABE  L  S ( 5 )  =  I TN2 
C 

C  CALCULATES  STRAIN  RATE  AND  TIME 

DO  10  J= 1 . NRR 
AE( J) =ESTRN1 (U+1  ) 

10  AT( J)=ELPST(U+1 ) *60 
GO  TO  3 

1  DO  11  J= 1 , NRR 
AE(U)=ESTRN2(U+ 1 ) 

11  AT( U) =ELPST(U+ 1 ) *60 
T  J  =  0 

3  DO  12  U=2 , NRR 

E(J)  =  (AE(J)-AE(J-1 ) )/(AT(U)-AT(J  -  1  ) ) 
TJ=(AL0G1O(AT(J))+AL0G1O(AT(J-1)))/2.O 

12  T(J)=10**TU 
I F ( MM  EO. 1 )  GO  TO  15 
XX=TWO 
GO  TO  16 

15  XX  =  ONE 

16  CONTINUE 
LABE  LS(9)=XX 

DETERMINE  MAXIMUM  AND  MINIMUM  VALUES 


U=  1 

25  d=J+1 

EMI N  =  E ( J ) 

EMAX  =  E ( J ) 

I F ( E ( U  )  .LE .0.0) GO  TO  25 
20  U=  J+  1 


^  V  i  HAi  1  W  £**}  i| 

*  ■  ■"»« 


■v  f  , 

«•' >  S  i 
W*f'  *U>^Jl*air* 


i  -v  fi  ,  i.-r 


«»  l  £T*ttr*rO  'tl*  * 


* 

K, 


78  1 

782 

783 

784 

785 

786 

787 

788 

789 

790 

79  1 

792 

793 

794 

795 

796 

797 

798 

799 

800 

80  1 

802 

803 

804 

805 

806 

807 

808 

809 

8  10 

8  1  1 

8  12 

8  13 

8  1  4 

8  15 

8  16 

8  17 

8  18 

8  19 

820 

82  1 

822 

823 

824 

825 

826 

827 

828 

829 

830 

83  1 

832 

833 

834 

835 

836 

837 

838 

839 

840 


o  n  o  non 


120 


IF(EMIN.LE  E(J)  OR  E(J)  LEO  0)GO  TO  21 
EM  I N=  E ( J ) 

21  I F ( EMAX  GE . E ( J) ) GO  TO  22 
EMAX  =  E  ( J  ) 

22  I F ( J  LT  NRR )  GO  TO  20 
DELETES  NEGATIVE  STRAIN  RATES 


<J  =  2 
1=2 

T  E ( I  ) =T ( J ) 

57  I  F  (  E  (  I  )  L  E  .0.0)  GO  TO  37 
1  =  1+1 
37  J  =  J+  1 

E  (  I  )  =  E  (  J  ) 

T  E (  I  )=T( J) 

IF ( J . LT . NRR )  GO  TO  57 
I  ND  =  0 
I ND= I  -  1 

IF(E(I).LE.O.)  IND= IND- 1 
I  F ( ANS2 .NE  . Y )  GO  TO  4 


PROMPT  FOR  STRAIN  RATE  -  TIME  PLOT  PARAMETERS 


202  FORMAT ( 

203  FORMAT ( 

204  FORMAT ( 

205  FORMAT( 


WRITE (6 . 202 ) 

WR I T  E ( 6 , 203 ) E  M I N 
WRITE(6.204)EMAX 
WR ITE(6,205)T( NRR ) 

FOLLOWING  ARE  THE  EXTREME  VALUES  IN  THE  T  VS  E  PLOT:') 
MINIMUM  STRAIN  RATE:  '.F12.6) 

MAXIMUM  STRAIN  RATE:  '.F12.6) 

FINAL  TIME :  ' , F 10. 3) 

WR I TE ( 6 , 206 ) 

206  FORMAT!'  INPUT  DESIRED  STRAIN  RATE  VALUE  AT  ORIGIN  OF  Y  AXIS 
•(TERMINATE  WITH  COMMA):') 

READ ( 5 . 302 )OPTNS(  1  1  ) 

302  FORMAT ( E  20 . 0 ) 

6  WR I T  E ( 6 , 207 ) 

207  FORMAT!'  INPUT  DESIRED  Y  AXIS  SCALE  ( UNI TS/ I NCH ) : ' ) 

READ! 5 , 303 ) OPTNS (  12) 

303  FORMAT ( E  20 . 0 ) 

IF (OPTNS! 12) .LT. lOOOOOO)  GO  TO  7 
WRITE (6 . 220) 

220  FORMAT! '•••  ERROR  ••*  ALL  INPUT  VALUES  MUST  BE  TERMINATED  BY  A 
•COMMA ' ) 


GO  TO  6 
7  WR I T E ( 6 , 208 ) 

208  FORMAT!'  INPUT  DESIRED  LENGTH  OF  Y  AXIS  (INCHES):') 

READ! 5 . 304 )OPTNS( 13 ) 

304  F  ORMAT ( E  20  O) 

WR I T  E ( 6 . 209 ) 

209  FORMAT!'  INPUT  DESIRED  X  AXIS  (TIME)  SCALE  ( UNI TS/ I NCH ) : ' ) 
READ ( 5 . 305 ) OPTNS ( 9 ) 

305  FORMAT ( E  20 . 0 ) 

WR I T  E ( 6 , 2 1 0 ) 

210  FORMAT!'  INPUT  DESIRED  LENGTH  OF  X  AXIS  (INCHES):') 

READ! 5 . 306 )OPTNS( 10) 

306  FORMAT ( E  20 . O ) 
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866 
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OP  T NS ( 6  )  =  1  O 
OP T NS ( 7  )  =  1  O 
OPTNS ( 16 ) =0 . 

OP T NS ( 2  1  )  =  1  .  O 
OPTNS ( 24 ) =BL ANK 
OPTNS ( 8 ) =0 
OPTNS (22) = 1 
OPTNS ( 23 ) =0 . 04 
LABELS(  13 ) =  BL ANK 
L  AB  E  L  S (  17 ) =BL  ANK 
C  ND=NRR-  1 

C 

CALL  STARTF  (GRAPH. TE ( 2  ). E ( 2  ).  I  NO , LABELS. -  1 . OPTNS ) 

C 

4  I F ( ANS 1  .  E  Q  .  Y  )  GO  TO  5 
GO  TO  60 

C 

C  CALCULATES  AXIS  VALUES  FOR  LOG  E  -  LOG  T  PLOT 
C 

5  XA= ALOG 10( T( 2 ) ) 

0PTNS(8)=FL0AT(  IFIX(XA)  )  -  1 
XE  =0 

XE=AL0G10(T( NRR ) ) 

XE=FLOAT ( IFIX(XE) )+ 1 
OPTNS ( 9 ) =XE-OPTNS(  8 ) 

XC* 10. 

YA=AL0G10(EMIN) 

OPTNS (  1  1  )= FLOAT (  I F I  X ( Y  A  )  )  -  1 
YE=AL0G10( EM  AX ) 

YE  =FLOAT ( IFIX(YE)  )+  1 
OPTNS (  12 ) = YE -OPTNS (  1  1  ) 

OPTNS ( 6 ) =2.0 
OPTNS ( 7 ) =  2 .0 
OPTNS (  1 7 )  =  1 .0 
OPTNS (  18) =2 .0 
0PTNS(2O)=O. 08 
OPTNS ( 2 1 ) = 1 .0 
OPTNS ( 2  2 )  =  1 
OPTNS ( 23 ) =0 . 04 
LABELS(  13 ) =LOG 
LABELS(  1  7 ) =LOG 
C 

CALL  STARTF  ( GR A PH . T E ( 2 )  . E ( 2  )  .  I ND . L AB E L S ,  1 . OP TNS ) 
CALCULATES  AND  PLOTS  BEST  FIT  LINES  FOR  LOG  E  -  LOG  T  PLOT 
K  =0 

40  K  =  K4- 1 

IF ( ANSM 1 . NE . Y )  GO  TO  41 
I F ( K . E 0 . 2 )  GO  TO  60 
NC  =  -  2 
GO  TO  42 

41  NC=K+ 1 
I  F ( K . EO  2 ) NC  =  - 2 

4  2  ND  =  4 

OPTNS ( 2  1  ) =2 .0 
OPTNS( 22  )  =4 
OPTNS ( 23 ) =0 . O 1 
D  =  0 

DX (  1  )  =  T ( 2  ) 


'  iv  Vt  ,'C j/  -AV 


>0  ■.  ;<:»  *«o 


<  •  t 4  .#  j  .  n  *£/  -  < .  >  *  y 


V- 


V 

• 

901 

902 

903 

904 

905 

906 

907 

908 

909 

9  10 

9  1  1 

9  12 

913 

9  1  4 

9  1  5 

9  16 

9  1  7 

9  18 

9  19 

920 

92  1 

922 

923 

924 

925 

926 

927 

928 

929 

930 

93  1 

932 
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D - BB 1 ( K . MM ) *  Al OG 10( DX (  1 ) ) +  BB0  ( K . MM ) 

DY(  1  )  =  10*  *D 
DX ( 2 ) =  T ( NRR ) 

D  =  BB 1 (K . MM ) » ALOG 10(DX ( 2  ) ) +BBO(K . MM ) 

D Y  (  2 ) = 10*  *D 
DY ( 3 ) =EMIN 

D=(AL0G1O(DY(3) )-BBO(K.MM) )/BB1 ( K . MM ) 
DX(3)=10*»D 
D Y ( 4 ) =  EMA  X 

D=(AL0G1O(DY(4)) -BBO(K, MM) )/BB 1 (K.MM) 

DX ( 4 )  =  10*  *D 
C 

DO  43  KK  =  1  . 4 

43  I F (DY ( KK ) . LT . O. OOOOOOI )  D Y ( KK ) =0  OOOOOO 1 


CALL  STARTF  ( GR A PH , DX , D Y . ND . L ABE L S . NC , OP TNS ) 
IF (K . LT . 2 )G0  TO  40 


60  I F ( MM . EO . 2 ) GO  TO  999 

IF ( XANS 1 . EO . Y )  GO  TO  50 
IF( IANS3 . EO. 1 )  GO  TO  999 
GO  TO  50 


999  CALL  STARTF  ( GR APH , TE ( 2 ) . E ( 2 ) . I ND . L AB E L S . O . OPTNS ) 


lOOO  RETURN 
END 


Sample  Run 


■  * 
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University  of  Alberta  -  Computing  Services  Device  V802  task  ^462 

*  s i g  gs22 
m  Password? 

? 

H  Term, Low, Internal/Teaching, Research 

*  Last  signon  was:  23:32:23 

*  User  ” GS2  2 "  signed  on  at  23:36:22  on  Wed  Ju 1  11/84 

«  $RUN  CPACK  4  =  ST  AR5  7  =  -P  T  =  4S 

THIS  PROGRAM  ACCEPTS  DATA  IN  THE  FOLLOWING  FORMAT: 

LABEL. TIME  CLOCK. 3  COUNT E R S , L 0 AD . C E L L  PR E S SUR E . 0 E F . 1 . D E F , 2 
IS  YOUR  DATA  IN  THIS  FORM?  (Y.N) 

NO 

DO  YOU  HAVE  AN  L  VDT  CALIBRATION  FACTOR? 

YES 

INPUT  THE  L  VDT  CALIBRATION  FACTOR  (MULTIPLIER)- 

O  02, 

INPUT  DI AME TER ( MM ) . LENGTH( MM ) . AND  WEIGHT(GM). 

SEPARATED  BY  COMMAS: 

54.76,49.2,0.0, 

INPUT  SAMPLE  NUMBER  (8  CHARACTER  MAX.): 

STARKE Y5 

INPUT  TEST  NUMBER  (8  CHARACTER  MAX.): 

STARKEY5 

INPUT  NUMBER  OF  ELEMENTS  IN  EACH  LINE  (+  COMMA)- 

10. 

INPUT  NUMBER  OF  STRAIN  GAUGES  (LVDTS):(1  OR  2): 

1  . 

THE  ONLY  DATA  THAT  THIS  PROGRAM  REQUIRES  FOR  INPUT 
IN  DEVICE  4,  IS  TIME  AND  DISPLACEMENT.  LIST  THE  ELEMENTS  IN 
ONE  LINE  OF  YOUR  INPUT  DATA,  IN  PROPER  ORDER.  USING  H R 1 ” 

FOR  THE  FIRST  STRAIN  GAUGE  READING,  "R2”  FOR  THE  SECOND  (IF 
THERE  IS  ONE).  AND  "TN" , "TD" . "TH" , "TM“ , ”TS" .  FOR  THE  TIME 
IN:  MONTHS. DAYS. HOURS. MINUTES. AND  SECONDS  (YOU  MAY  HAVE  ONLY 
ONE  OF  THESE  TIME  PARAMETERS).  USE  THE  DUMMY  VARIABLE  "DU" 

FOR  ALL  OTHER  VALUES.  (SEPARATE  THE  ELEMENTS  WITH  COMMAS) 

DU . DU . DU . DU . TH . DU . R 1 . DU . DU , DU . 

WOULD  YOU  LIKE  THE  DECELERATING  CREEP  DATA 
PLOTTED  BEFORE  THE  PROGRAM  ATTEMPTS  TO  ISOLATE 
ACCELERATING  CREEP? 

NO 

WOULD  YOU  LIKE  THE  PROGRAM  TO  EXAMINE  THE 
ACCELERATING  DATA  FOR  AN  OPTIMUM  FIT? 

YES 

INPUT  DESIRED  LIMIT  FOR  TEST  OF  SLOPE  SIGNIFICANCE  IF 
OTHER  THAN  10.  (REAL  NUMBER.  TERMINATED  WITH  A  COMMA): 

WOULD  YOU  LIKE  A  LOG  STRAIN  RATE  -  LOG  TIME  PLOT?(Y,N) 

YES 

WOULD  YOU  LIKE  A  STRAIN  RATE  -  TIME  PLOT?(Y.N) 

YES 

FOLLOWING  ARE  THE  EXTREME  VALUES  IN  THE  T  VS  E  PLOT: 

MINIMUM  STRAIN  RATE:  0.028226 

MAXIMUM  STRAIN  RATE:  3.387532 

FINAL  TIME:  228028  500 

INPUT  DESIRED  STRAIN  RATE  VALUE  AT  ORIGIN  OF  Y  AXIS  (TERMINATE  WITH  COMMA) 

0.0. 

INPUT  DESIRED  Y  AXIS  SCALE  (UNITS/INCH) 

0.4  . 

INPUT  DESIRED  LENGTH  OF  Y  AXIS  (INCHES): 

10.  . 
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INPUT  DESIRED  X  AXIS  (TIME)  SCALE  (UNITS/INCH) 
25000  . 

INPUT  DESIRED  LENGTH  OF  X  AXIS  (INCHES) 

10.  . 


GRAPH  preliminary  VERSION  DEC.  1.  1982 

PLOT  FILE  NAME  IS  -PDF 


SUMMARY  file  name  is  -summary 


T 

T 


—  T 
TT  + 
TT| 

—  T 

=  *f 

|“l  ' 

|  =  =  T 

♦♦K 

T  +  TT 


OPTION 

1.  PLOT 

2.  BLOW-UP 

3.  REDRAW 

4.  SUBPICTURES 

5.  MTS-SDS 

6.  CONTINUE 


T-  ! 

--T=_T_-_f  _  +  _  -  -  =- 

=  =  + _ +T_=  =  =  +  =  +  =_  + _ =  =  -t-TT  =  T  =  T  =  +  +  +  + 

T-  -T-'  ->T-«  ->T->  -T++TT- - T  -  -T-*  ->T  ->T~  ^T~* 

I+  +  T+T  +TTT  +  T  |  TTT  +  T-’T  + 


SELECT  MENU  OPTION 

?  1 


SELECT  MENU  OPTION 

?  6 


TTT  + 
T  =  + 

T  =  + 
TT  T 
*+  T 
*T  = 


T 

+  +  +  -t- 


+  •*■  |  -i  + 

■_  -  TTT  -_«■*-  +T 

->  I  T  =  =  T  - 

|+||-  T+++ 

”,T  |  T  +|| 

- _ T  |  T  =  — 

+  T  I 
I  + 

*  +  +  4  +  + 


1  . 
2  . 

3  . 

4  . 

5. 

6. 


OPTION 

PLOT 

BLOW-UP 

REDRAW 

SUBPICTURES 

MTS-SDS 

CONTINUE 
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T  |  TT 


I-  T  T  ->  T  T  T  T  T  T  =  |  T-  T  T  ->  T  T  T 
|  +  +  T ♦ T  +TTT*T|TTT  +  T-'T* 


SELECT  MENU  OPTION 

?  1 


SELECT  MENU  OPT  I  ON 

?  6 
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Input  File 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  1 

12 

13 

14 

15 

16 

1  7 

18 

19 

20 

2  1 

22 

23 

24 

25 

26 

27 

28 

29 

30 

3  1 

32 

33 

34 

35 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

46 

47 

48 

49 

50 

5  1 

52 

53 

54 

55 

56 

57 

58 

59 

60 
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0.0. 0.0. 12. .0.969  .0. 1241 1 . .0. 
0.0. 0.0. 16. .0.982.5,0. 13790. .0. 
0.0.0,0.33. .0.988, .0. 13445.25.0. 
0.0.0.0.39  .0.991 . ,0. 13790. .0. 
0.0.0,0.58. .0.992.8.0. 13790. .0. 

0. 0.0.0. 64 . .0, 993 . 7 .0. 13790  .0. 
0.0.0.0,66. .0.994, .0. 13790. .O. 
0.0.0.0.71 . .0.994,5.0. 13790. .0. 
0,0.0.0.82. .0.995. .0. 13790. .0. 
0.0.0.0.89.  .0.996.  .0.  13790.  .0. 
0.0.0.0.95. .0,996 .2.0, 13790. .0. 
0.0. 0.0. 106  .0,997.2.0. 13790. .0. 
0.0. 0.0. 1 12 . .0, 998  3.0.13790  .0. 
0.0. 0.0. 129  .0,9988.0. 13790. .0. 
0.0. 0.0. 136  .0, 1 000  . 0 . 13790  . 0 . 
0.0. 0.0. 182  .0. 1002 .3.0. 13790  .0. 
0.0. 0.0. 206  .0. 1003  2.0. 13790. .0. 
0.0.0.0.21 1 . .0. 1004  7,0. 13790. .0. 
0.0.0.0,226  .0. 1005 . 1.0. 13790  .0. 
0.0.0,0.235. .0. 1005.5.0, 13790  .0, 
0,0,0.0,274. .0. 1007  2.0. 13790. .0. 
0.0.0.0.281 . .0. 1008 .8.0. 13790  .0. 
0.0.0.0.326, ,0. 1009  8.0. 13790  .0. 
0.0.0.0.346. .0. 1011. .0. 13790. .0. 
0.0.0.0.370. .0. 1012. .0. 13790. .0. 
0.0.0.0.382. .0. 1012.2,0, 13790. .0. 
0.0.0.0.393. .0. 1013. ,0. 13790  .0. 
0.0.0.0.401 . ,0. 10138.0. 13790  .0. 
0,0,0,0,417. .0. 1014. .0. 13790. .0. 
0.0.0.0,522. .0. 1017  .0. 13790. .0. 
0.0.0.0.537. ,0. 1018. ,0. 13790. .0. 
0.0.0,0.561 . .0. 1018.4.0. 13790. .0. 
0.0.0,0,585. .0. 1019. .0. 13790. .0. 
0.0,0.0,591. .0. 1019. .0. 13790. .0, 
0.0.0,0.610. .0. 1019.8.0. 13790. ,0. 
0,0.0.0.642. .0. 1020.5,0. 13790. .0. 
0.0.0.0.686. .0. 1022. .0. 13790. .0. 

0 . 0 ,0,0.706 . .O. 1022.5.0. 13790. .0. 
0.0.0.0.736. ,0, 1023. 1.0. 13790  .0. 
0.0.0.0.808 . .0. 1025 . ,0. 13790. . 0 . 
0.0.0.0.849, ,0. 1026.2,0. 13790. .0. 
0.0.0.0.876. .0. 1027,2.0. 13790. .0. 
0.0.0.0.973. .0. 1028.8.0. 13790. .0. 
0,0. 0.0. 1020  .0. 1029. 1 ,0. 13790. .0. 
0.0. 0.0. 1071 . .0. 1030. 9.0. 13790. .0, 
0.0. 0.0. 1185. .0. 1032  2.0. 13790. .0. 
0.0. 0.0, 1210. .0. 1032.9,0. 13790. .0. 
0.0. 0.0. 1264. .0. 1033.2.0. 13790. .0, 
0.0. 0.0. 1281. .0. 1033.6,0. 13790. .0. 
0.0. 0.0. 1353. .0. 1033.4.0. 13790. .0. 
0.0. 0.0. 1381. ,0. 1033.3.0. 13790. .0. 
0.0. 0.0. 1402 . .0. 1034 , 1.0. 13790. ,0. 
0.0. 0.0.  1449.  .0.  1035.9.0,  13790.  .0. 
0.0. 0.0.  1523.  .0.  1037,6.0.  13790.  .0. 
0.0. 0.0. 1545. .0. 1037.6.0. 13790  .0. 
0.0. 0.0. 1571 . .0. 1037  8.0. 13790. .0. 
0.0. 0.0. 1593  .0. 1038.2.0. 13790. .0. 
0.0. 0.0.  1617.  .0.  10386.0.  137  90.  .0. 
0.0. 0.0.  1693.  .0.  1039  6.0,  137  90  .0. 
0.0. 0.0.  1762.  .0,  1039.9,0.  137  90  .0. 


). 

u  Of  -  c  .  J  U.o 
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6  1 

62 

63 

6  4 

65 

66 

67 

68 

69 

70 

7  t 

72 

73 

74 

75 

76 

77 

78 

79 

80 

8  1 

82 

83 

84 

85 

86 

87 

88 

89 

90 

9  1 

92 

93 

94 

95 

96 

97 


129 


0.0.0.0.1785  .0.1040  1.0.13790  .0 
0.0. 0.0. 1881 . .0. 10409.0. 13790. .0, 
0.0. 0.0. 1905. .0. 1041  2.0. 13790. .0, 
0.0. 0.0. 1929  .0. 1041  7.0. 13790. .0. 
0.0. 0.0. 1953  .O. 1042  .0. 13790. .0. 
0.0.0.0.2025. .0. 1042  9.0. 13790. .0. 
0. 0. 0. 0. 2049 , .0. 1043  1.0. 13790 . . 0. 
0.0.0.0.2073  .0. 1043  1.0, 13790. .0. 
0,0,0.0,2121. .0. 1043. 1.0. 13790. .0. 
0.0.0.0.2193.  .0.  1042.7.0.  137  90.  .0.' 
0.0.0.0.2216  .0. 104 26.0. 13790. .0. 
0.0,0.0.2240. .O. 1043  .0. 13790  .0. 
0.0.0.0.2265  .0. 1042.8.0. 13790. .0. 
0.0.0.0,2361. .0. 1043  2.0. 13790. .0, 
0.0.0.0.2385  .0. 1043.3.0. 13790  .0. 
0.0,0.0.2409. .0. 1043. .0. 13790. .0. 
0.0.0.0.2439. ,0. 1042.7,0, 13790. ,0. 
0.0.0.0.2460. .0. 1042.6.0. 13790. .0. 
0.0.0.0.2540. ,0. 1043.8.0, 13790. ,0. 
0.0.0,0.2552. .0. 1044. .0, 13790. .0. 
0.0.0.0.2580. ,0, 1044.2.0. 13790. .0. 
0,0.0.0,2601 . .0. 1044.3.0, 13790. .0. 
0,0.0.0,2722.5.0. 1046. .0. 13790. .0. 
0,0.0.0.2890. .0. 1047. ,0. 13790. .0. 
0. 0 , 0. 0 . 3009  .0. 1050 .6.0. 13790. .0. 
0 , 0. 0 . 0. 3056 .  . 0  .  105  2 .  .0.  137  90.  . 0 . 
0.0,0.0.3082 . .0. 1052.8.0. 13790 . .0. 
0.0.0,0.3201  .0. 1058.7,0. 13790. .0. 
0.0.0.0.3272. .0. 1066. .0. 13790. ,0. 
0.0,0.0,3371. .0. 1072.4.0. 13790. .0. 
0.0.0.0.3466 . .0. 1073 . .0. 13790. .0. 
0.0.0.0.3560. .0. 1075.2.0. 13790. .0. 
0.0.0.0.3610. .0. 1082.4.0. 13790. ,0, 
0.0.0.0.3661 . .0, 1084.5,0, 13790. ,0. 
0.0.0.0.3730. .0. 1088. .0. 13790  .0. 
0.0.0.0.3800  ,0.  1090.  .0.  13790.  .0. 
0.0,0.0,3825. .0. 1092.7,0. 13790. . 1. 


Output  File 
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TEST  NUMEER  -  STARKEY5 

SAMPLE  NUMBER  =  STARKEY5 

SAMPLE  LENGTH  =  49.200 

SAMPLE  DIAMETER  =  54  760 

SAMPLE  WEIGHT  =  0.0 

INITIAL  LOAD 
INITIAL  CELL  PRESSURE 
INITIAL  READING  FOR  LVDT1  = 
INITIAL  READING  FOR  L  VD  T  2  = 


NUMBER 

TIME 

LOAD 

HOURS 

KN 

1 

0.0 

0 . 0 

2 

4.000 

0.0 

3 

2 1 . OOO 

0 . 0 

4 

27 .OOO 

0.0 

5 

4  6 . OOO 

0.0 

6 

52.000 

0.0 

7 

54 . OOO 

0.0 

8 

59 . OOO 

0.0 

9 

70 . OOO 

0.0 

10 

77.000 

0.0 

1  1 

83.000 

0.0 

12 

94  OOO 

0.0 

13 

1 00 . OOO 

0.0 

14 

117. OOO 

0.0 

15 

124 .000 

0.0 

16 

170.000 

0.0 

1  7 

194 .OOO 

0 . 0 

18 

199 . OOO 

0.0 

19 

214. OOO 

0.0 

20 

223  OOO 

0.0 

2  1 

262 .OOO 

0.0 

22 

269.000 

0.0 

23 

314. OOO 

0.0 

24 

334  OOO 

0.0 

25 

358 . OOO 

0.0 

26 

370.000 

0.0 

27 

38 1 . OOO 

0.0 

28 

389 . OOO 

0.0 

29 

405 . OOO 

0  0 

30 

510. OOO 

0.0 

3  1 

525 . OOO 

0.0 

32 

549  OOO 

0 . 0 

33 

573. OOO 

0.0 

34 

579 , OOO 

0.0 

35 

598 .000 

0  0 

36 

630.000 

0.0 

37 

674 . OOO 

0.0 

38 

694 . OOO 

0.0 

39 

724 . OOO 

0.0 

40 

796 . OOO 

0.0 

4  1 

837 . OOO 

0.0 

42 

864  OOO 

0.0 

43 

96 1 . OOO 

0.0 

44 

1008  OOO 

0.0 

mm 

mm 

gm 

0.0  KN 
0  0  KP  A 
969.0000000000  mm 
O  .  O  mm 

CELL  PRESSURE  ENGG  STRAIN 
KPA  (MICRO) 


0.0 

0.0 

0.0 

5487.801 

0.0 

7723 . 574 

0.0 

8943 .086 

0.0 

9674 .813 

0.0 

10040.629 

0.0 

10162. 598 

0.0 

10365  852 

0.0 

10569 . 105 

0.0 

10975 . 609 

0.0 

1  1056 . 887 

0.0 

1  1463 . 395 

0.0 

1 1910.586 

0.0 

12113  836 

0.0 

12601 . 62 1 

0.0 

13536  602 

0.0 

13902 . 4 18 

0.0 

145  12.  172 

0.0 

14674 .832 

0.0 

14837 . 395 

0  0 

15528  434 

0.0 

16178. 879 

0.0 

16585 . 383 

0.0 

17073 . 168 

0.0 

17479 . 672 

0.0 

17560. 953 

0.0 

17886  .  176 

0.0 

18211. 398 

0.0 

18292 . 680 

0.0 

19512 . 191 

0.0 

19918. 695 

0.0 

2008 1 . 250 

0.0 

20325 . 199 

0.0 

20325 . 199 

0.0 

20650 . 422 

0.0 

20934 . 94 1 

0.0 

2 1544 . 699 

0.0 

21747.938 

0.0 

21991. 898 

0  0 

22764  219 

0.0 

2325  1 . 996 

0  0 

23658 .496 

0.0 

24308 . 938 

0.0 

24430.918 

ENGG.  STRAIN 
*2  (MICRO) 
0.0 
0  0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
O  0 

o .  o 
0.0 
0  o 

O  0 
0.0 
0.0 
0.0 
0.0 
o  o 
0.0 
0.0 
0  0 
0.0 
0.0 
O  0 
0.0 
0.0 


' 


ii  & 


. 

. 


t 

/ 

c 

- 
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45 

1059 . OOO 

46 

1173. 000 

4  7 

1 198 .OOO 

48 

1252 . OOO 

49 

1269 . OOO 

50 

1341. OOO 

5  1 

1369 .000 

52 

1390. OOO 

53 

1437 . OOO 

54 

1511. OOO 

55 

1533.000 

56 

1559  OOO 

57 

1581. OOO 

58 

1605 . OOO 

59 

1681. OOO 

60 

1750.000 

6  1 

1773  OOO 

62 

1869 . OOO 

63 

1893  OOO 

64 

1917. OOO 

65 

1941. OOO 

66 

2013. OOO 

67 

2037 . OOO 

68 

206 1 . OOO 

69 

2109.000 

70 

2181. OOO 

7  1 

2204 . OOO 

72 

2228  OOO 

73 

2253 . OOO 

74 

2349 . OOO 

75 

2373 . OOO 

76 

2397 . OOO 

77 

2427 . OOO 

78 

2448 .000 

79 

2528 . OOO 

80 

2540.000 

8  1 

2568 . OOO 

82 

2589 . OOO 

83 

27 10.500 

84 

2878 . OOO 

85 

2997 . OOO 

86 

3044 . OOO 

87 

3070 . OOO 

88 

3189  OOO 

89 

3260 . OOO 

90 

3359 . OOO 

9  1 

3454 . OOO 

92 

3548  OOO 

93 

3598 . OOO 

94 

3649  OOO 

95 

3718  OOO 

96 

3788  OOO 

97 

3813  OOO 

0  0 

0.0 

0.0 

0  0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0 . 0 

0.0 

0  0 

0  0 

0.0 

0.0 

0 . 0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0  0 

0.0 

0  0 

0  0 

0.0 

0.0 

0  0 

0.0 

0.0 

0  0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0  0 

0.0 

0.0 

0.0 

0  0 

0.0 

0  0 

0.0 

0  0 

0  0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0  0 

FIT  OF  DECELERATING  CREEP  DATA 


25162  559 

0  0 

2569 1016 

0  0 

25975 . 559 

0  0 

26097  516 

0  0 

26260  176 

0.0 

26178  797 

0  0 

26138  219 

0.0 

26463  438 

0  0 

27195  059 

0.0 

27886  199 

0.0 

27886 . 199 

0.0 

27967  477 

0.0 

28 1 30  039 

0.0 

28292  699 

0  0 

28699  219 

0 . 0 

28821  078 

0  0 

28902  457 

0  0 

29227  598 

0  0 

29349 . 559 

0 . 0 

29552  816 

0.0 

29674  777 

0.0 

30040.598  ■* 

0.0 

30121  957 

0  0 

30121  957 

0.0 

30121 . 957 

0.0 

29959  297 

0.0 

29918 .719 

0.0 

3008 1.277 

0.0 

29999 . 996 

0.0 

30162  578 

0  0 

30203 . 258 

0.0 

3008 1.277 

0 . 0 

29959 . 297 

0.0 

29918  719 

0.0 

30406 . 496 

0.0 

30487  797 

0.0 

30569 . 059 

0 . 0 

30609 . 758 

0.0 

3 1 300 .797 

0  0 

31707 . 297 

0  0 

33170.738 

0.0 

33739  816 

0.0 

34065 . 039 

0.0 

36463  379 

0 . 0 

39430 . 879 

0.0 

42032 . 457 

0.0 

42276 .4  18 

0.0 

43170.695 

0  0 

46097  496 

0.0 

4695  1.195 

0.0 

48373 . 977 

0.0 

49186  977 

0.0 

50284 . 496 

0.0 

TO  POWER  LAW  ****** 

DATA  FROM  LVDT  NO.  1 
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TRANSFORMED  DATA 


T  I  ME 

STR  RATE.  E 

(MIN) 

(MICRO  E/MIN) 

5J9  9077 

2  .  19  1935 

1428  7019 

3  387532 

2114.5154 

0  64  1865 

2934  4822 

1 .016156 

3179  4292 

1 . 016406 

3386 . 6707 

0.677513 

3855  8955 

0 . 307960 

4404  9922 

0 . 967866 

4796  6172 

0 . 225770 

5299  7188 

0  6 1592  1 

5817  2070 

1  242198 

6489  9844 

0. 199265 

7226  9219 

1.161 393 

8711.3477 

0 . 33876  1 

10896 .2109 

0.254039 

1 1789 .0313 

2 .0325  13 

1238  1  . 8086 

0. 180734 

13 107 . 1992 

0. 301042 

14502 . 8438 

0. 295316 

15928 . 5742 

1 . 548678 

17437 . 7852 

0. 150557 

19430  6875 

0 . 406488 

20747 .4570 

0. 282294 

2 1836 . 9688 

0.  1 1289  1 

22527 . 5273 

0 . 492762 

23098 . 6875 

0  677547 

238  15.  1289 

0 . 084668 

27268 . 6094 

0. 193573 

31046. 6602 

0.451671 

322 11.9141 

0.  1  12885 

33740 .1914 

0. 150586 

352  13. 703  1 

0 . 24638  1 

36827 .3516 

0. 148187 

39097 . 6328 

0. 230969 

4 1035 .5195 

0.  169365 

42530 .3516 

0.  135534 

45548  6719 

0.  178778 

48974 , 5039 

0.  198283 

51023  4766 

0.250926 

54672 . 4805 

0  111  760 

59053 .0117 

0.043255  • 

61990.9844 

0.239098  ■ 

66872 . 3750 

0.077260  • 

71125. 8750 

0  189695  - 

73482 . 0000 

0.037641  - 

773 1 7 . 6875 

0.021273  - 

81467 . 4375 

0.070954  - 

84798 .  1875 

0.259440  - 

88733 .3125 

0 .  1355  18  - 

92423 . 2500 

0.036611  - 

94197. 3750 

0.123153  - 

95576 . 9375 

0 . 1 12958  - 

98553  3125 

0  089149  - 

102908 .8750 

0.029435  - 

105687 . 5000 

0.058970  - 

LOG  E 

O  34  08  28  E  +  00 
0. 529883E + 00 
-0,  1 9  2  5  5  6  E  *  00 
0  696047E -02 
O.  706724E -02 
-0.  1 6908  2  E  +00 
-0  5  1  1505E  +00 
-0  14  1846E-01 
-O . 646333E+  00 
-0 . 2  10475E+00 
0 . 94 1908E -01 
-0 . 700570E+00 
O . 649793E -0 1 
-0.470107E+00 
-O. 595099E+00 
0 . 308033E+00 
-0.74296 1 E  +00 
-0.52 1373E+00 
-O. 5297 13E+00 
O  .  18996 1  E  +  00 
-0.822299E+00 
-0 . 390953E+00 
-0. 549298E  +00 
-0 . 94  7 3 4  2  E  +00 
-0 . 307363E+0O 
-0. 169060E+00 
-O. 107228E+01 
-O. 7 13 155E+00 
-0  3451 78 E+00 
-0 . 947363E  +00 
-0.822215E+00 
-0 . 608  39  3E  +00 
-O. 829 189E+00 
-0 . 636447E+- 00 
-0.77  11 76  E  +00 
-O. 867952E+00 
-O.  747686E+-00 
-0. 7027 13E+00 
-O. 600454E+00 
-0.95  17 1 5  E  +00 
-0.  1  36396E+0 1 
-0 . 62 1424E+00 
-0.11  1205E+-01 
-0.72  1943E+-00 
-0.  1 424  34  E  +0 1 
-O. 1672 16E+01 
-0. 1 14902E+01 
•O  .  5 8 596 3 E  +  00 
-O  .  868004 E  +-00 
-0.  1 4  36  38  E  +0 1 
O. 909553E+00 
0.947082E+00  - 
0.  104988E  +0 1  - 
0.  1531  14E+01  - 
O.  1 2  29  37  E  +0 1  - 


LOG  EE 


LOG  E  -  LOG  EE 


O.  856889E+00 
O. 44 1976E+00 
0.27 1600E  +00 
0  129 189E+00 
O. 943499E-01 
O  669088E-01 
O .  1 05209  E -O 1 
-O . 473347E -0 1 
-O. 843477E-01 
-0.  1 2  7  6  9  3  E  +  00 
-0 . 168 180E+00 
-O. 2 15739E+00 
-O. 262479E  +  00 
-0 . 34366 1 E+00 
-O  .4409  1  2E+-00 
-0.4751 37E+00 
-0 . 496456E+-00 
-0.521  197E+-00 
-0. 565 168E+00 
-O . 6059 18E+00 
-0.645257E+00 
-0 . 692284E-+QO 
-O.  720778E+-00 
-0. 743019E+-00 
-O . 756549E  +  00 
-0. 767429E+00 
-O.  780704  E-+00 
-0 .83955  1  E+-00 
-O . 895938E  +00 
-091 1949E+00 
-0 . 932094  E  +  00 
-O. 950669E+00 
-0. 970140E+00 
-O  996137E+00 
-O . 10 1 7 16E+01 
-O. 10327 1E+01 
-O.  106250E  +0 1 
-0.  1 09402  E  +0 1 
-0.111 183E+01 
-0. 114 185E+01 
-O  1  17534E+01 
-O.  1  1 964  4  E  +0 1 
-O.  1 22938E  +0 1 
-O. 1256 18E+01 
-O.  1 2  7034  E  +0 1 
-O.  129245E  +  01 
-O. 13 15 17E+01 
-0 .  1 33258  E  +0 1 
-O. 135230E+01 
-O.  1 370OOE  +0 1 
-O.  1 37826E  +  0 1 
O.  1  38458E  +0 1 
O.  1  3979 1 E  +  0 1 
0.14  1670E+01  - 
O. 142828E+0 1 


-0  51606 1 E+00 
O  8  7  907 7 E -0 1 
-0  464 156E+00 
-O.  1  222  29E  +00 
-O . 872826E -O 1 
-O . 23599 1 E+00 
-0 . 522026E  +00 
0  33 1501E-01 
-0.56 1985E+00 
-O . 8278 19E -0 1 
O.  262370E+00 
-0. 484830E+00 
O  327458E  +00 
-O.  126445E  +  00 
-0 .  154  1 87  E  +  00 
O  783170E+00 
-0. 246505E+00 
-O.  175893E-03 
O.  354553E-01 
O.  795879E+00 
-O.  1 7  704  2  E  +00 
O. 301331 E+00 
O. 17 1480E+OO 
-O. 204323E+00 
O. 449 186E+00 
0. 598369E  +  00 
-0. 291577E+00 
0.  1 2  6  3  9  6  E  +  00 
O . 550760E  +00 
-0. 354 140E-01 
0. 1098 78 E+00 
0 . 342276E  +  00 
0. 14095 1E+00 
0 . 359690E  +  00 
0  245984E+00 
0. 164756E+00 
0.314817E+00 
0. 39 1 304  E  +  00 
0.51 1375E+00 
O.  1 90 1 32  E  +00 
-O  1886 18E+00 
0 . 57  50 1 7 E  +  00 
O. 1 1 7332E+00 
O . 534  2  34  E  +00 
-O.  1 5  3  9  9  8  E  +  00 
-O . 3797  1  4  E  +  00 
0. 166 145E+00 
0. 74662 lE+OO 
O. 484292E+00 
-0.663815E-01 
O.  4687  1  1 E  +  00 
0 . 4  3  7  502  E  +00 
O . 34  8026  E  +00 
0  1 1 4  4  38  E  +00 
O  1989 14E+00 


W 

1 .00 
1 .00 
1  .OO 
1  OO 
1 .00 
1 . 00 
1  OO 
1 . 00 
1 . 00 
1 .00 
1 . 00 
1  .00 
1  .  OO 
1  .00 
1  OO 
1  .OO 
1  .OO 
1  .OO 
1 .00 
1  OO 
1 . 00 
1  .00 
1  .OO 
1 . 00 
1 . 00 
1  .OO 
1 .00 
1  .OO 
1 .00 
1 . 00 
1  .  50 
1  .25 
1.13 
1 . 06 
1 .03 
1 . 02 
1.01 
1 .00 
1  OO 
1 . 00 
1 . 00 
1 . 00 
1 .00 
1  .OO 
1  .00 
2.00 
1 . 50 
1 . 25 
1 . 63 
1.31 
1  .  16 
1 .08 
1 . 04 
1  02 
1.01 
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10922 1  7500 
112857  4375 
1 14297  5000 
1  15737 .3  125 
1 18599  7500 
1249 17.1 250 
135628  0625 


0  056448 
0 . 08  1695 
0  141151 
0 . 084695 
0  08468 1 
0  000692 
0.007782 


-0. 124835E*01 
-0  1072 1 4  E  4  0 1 
-0  8  503 1 5E  4 00 
-0.  1072 1 4  E  4  0 1 
-0. 107222E+0 1 
-0. 3 15963E+01 
-0. 2 1 0889E  40 1 


-0. 14425BE401 
-0  14  568 1 E  4 0 1 
*0  146232E*01 
-0.  146776E401 
-0.  147837E+01 
-0 . 15OO93E+01 
-0.  1  53668E  40 1 


0.  1  9  4  2  2  4  E  4  OO 
O.  384664E+  OO 
0  6 1 2000E  4  00 
O. 3956 1 4  E  4  00 
O  406 1 57E  400 
-O.  16587 1 E  4  0 1 
-0. 5722 1 3E  4 00 


1  00 
1  .00 
1  OO 
1  00 
1  .00 
4  00 
3 . 50 


FIT  PARAMETERS 
INTERCEPT.  BO 
SLOPE.  B 1 

CONFIDENCE  LIMIT  ON  BO  T. 

CONFIDENCE  LIMIT  ON  B1  T 

DURBIN  WATSON  STATISTIC 
DEGREE  OF  FREEDOM  FOR  DW 
UPPER  LIMIT  FOR  DW 
TEST  OF  SLOPE  SIGNIFICANCE 


O  359890E  +0 1 
-  1  0006285 
0  5185753 
0  1145718 
1  677 
75 

1  650 
76 . 249 


DATA  FOR  COMPARISON  TESTS 
WEIGHTING  69.5 

MEAN  STRAIN  -0.895  MEAN  TIME  4.491 

SSDY  43.673  SSDX  22.827 

SPDXY  -22.835  SSDYX  20.825 

CHECK 

SUM  OF  RESIDUALS  -0.001285 


DATA  STARTS  AT  *  1 

DATA  ENDS  AT  *  75 

FIT  PARAMETERS  FROM  PREVIOUS  ITERATION 

TEST  OF  SLOPE  SIGNIFICANCE  77.291 

DURBIN  WATSON  STATISTIC  1.500 


FIT  OF  ACCELERATING  CREEP  DATA  TO  POWER  LAW 


PARAMETERS  FOR  EVALUATING  TOTAL  FIT.  LVDT  NO.:  1 


RANGE 

ACCELERATING  creep 
INTERCEPT, BO  SLOPE. B1 

R 

DW 

TOTAL 

R 

F  I  T 

56  - 

96 

-031 1423E+02 

5  7150812 

44  2  10 

1  .  588 

46 

49  1 

57  - 

96 

-0  349556E  +02 

6  4425755 

62  142 

1  .  844 

47 

.  474 

60  - 

96 

-0 . 4 1 7652E+02 

7 . 73994  16 

90 . 227 

1.617 

48 

796 

6  1  - 

96 

-0 . 448557E402 

8  3275833 

112.218 

2  630 

49 

2  15 

62  - 

96 

-0. 492026E+02 

9 .  15338  13 

197 . 507 

1  609 

49 

609 

63  - 

96 

-0 . 440275E+02 

8 . 1738224 

57  .  144 

1.681 

49 

.  039 

64  - 

96 

-0 . 435989E+02 

8 . 0928059 

52.815 

1  .  690 

48 

973 

DATA  FROM  LVDT  NO.  1 


TRANSFORMED  DATA 


' 
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time 

(MIN) 

149530 .1875 
167579 .1250 
176213  5000 
18  1223 . 9375 
183417.8125 
187735 . 7500 
193457 . 9375 
198547 .  1875 
204369  5000 
2 10040. 2500 
2  14373 .8  125 
217404. 0000 
220999  5000 
225169  7500 
228028 . 5000 


STR  RATE  .  E 
(MICRO . E/MIN) 

O  019336 
0.016897 
O.  1825  1  1 
O.  179970 
0 . 187042 
O.  3  14885 
0.676  174 
0.4 18027 
0.023469 
O. 139765 
O . 957  12  1 
0. 260870 
0  325762 
0. 176027 
0.714471 


LOG  E 

-0.17 1364E+01 
-O. 177220E+01 
-O  7387  1  lE+OO 
-O. 744800E+00 
-O. 72806 1E+00 
-O. 50 1 848E  *00 
-0. 169942E+00 
-O . 3  78  7  96  E  *  00 
-O.  1 62950E  +  0 1 
-0. 854600E  +00 
-O. 190331E-01 
-0 . 583576E  +00 
-O. 487099E+00 
-0. 75442 1 E  +00 
-0. 146015E+00 


LOG  EE  L 

-O.  172083E+01 
-O. 13203 1E+01 
-O  1  14372E+01 
-O  1045 18E+01 
-0  100290E+01 
-0.9211 12E+00 
-0  8  15582E+00 
-O  7  2  4  3  1 9 E  +  00 
-O. 622742E+00 
-0. 526550E+00 
-0  4 54 7 7 3 E  +00 
-O . 40544 1 E+00 
-0  34777BE+00 
-O  28  207  4  E  *00 
-0. 237747E+  00 


LOG  E  -  LOG  EE 

1  O  718498E-02 
1  -0. 45 1883E+00 
1  O  405012E+00 
1  O . 30038 1 E+00 
1  0.274838E  +  00 

0  O  419264E+00 
0  0.64564 1E  +  00 

0  0. 345523E  *00 

0  -O.  1 006  7  6  E  *  O 1 
0  -0 . 328050E+00 
9  0.435740E+00 

9  -O  17  8  1  35  E  +00 
9  -0  13932 1 E  +00 
9  -0.472347E+00 
9  0.917320E-01 


W 

9 . 50 
5 . 25 

3.13 
2 . 06 
1 . 53 
1  .  27 

1.13 
1  07 
1 . 03 
1  .02 
1.01 
1 .00 
1  00 
1  00 
1 . 00 


FIT  parameters 

INTERCEPT.  BO 
SLOPE.  B  1 

CONFIDENCE  LIMIT  ON  BO  T 

CONFIDENCE  LIMIT  ON  B1  T. 

DURBIN  WATSON  STATISTIC 
DEGREE  OF  FREEDOM  FOR  DW 
UPPER  LIMIT  FOR  DW 
TEST  OF  SLOPE  SIGNIFICANCE 


-0  435989E+02 
8.0928059 
5 . 86 10058 
1 . 1 169233 
1 .690 
33 

1  .  505 
52.815 


DATA  FOR  COMPARISON  TESTS 
WEIGHTING  30.0 

MEAN  STRAIN  -1.135  MEAN  TIME  5.247 

SSDY  11.981  SSDX  0116 

SPDXY  0.943  SSDYX  4  340 

CHECK 

SUM  OF  RESIDUALS  0.001468 


DATA  STARTS  AT  M  64 
DATA  ENDS  AT  #  96 

FIT  PARAMETERS  FROM  PREVIOUS  ITERATION 

TEST  OF  SLOPE  SIGNIFICANCE  74.130 

DURBIN  WATSON  STATISTIC  1  412 
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ADDITIONAL  TEST  RESULTS 
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LOG  STRRIN  ROTE 
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log  Tire 


Axial  Stress  =  6.38  MPa 


Figure  B.1  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  U9~9#11 
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Axial  Stress  =  6.38  MPa 


Figure  B.2  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  U9~9#21 
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Axial  Stress  =  7.18  MPa 


Figure  B.3  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  U10-9#2 


p  - 
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Figure  B.4  Logarithm  plot  of  strain  rate  (micro-st ra in/min ) 
vs  time  (min)  Test  S9-9#7 
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Axial  Stress  =  3.99  MPa 

Figure  B.5  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  S9-9#15 
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Axial  Stress  =  3.99  MPa 


Figure  B.6  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  S9-9#17 
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Axial  Stress  =  3.99  MPa 


Figure  B.7  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  S10-9#2 
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Axial  Stress  =  2.39  MPa 


Figure  B.8  Logarithm  plot  of  strain  rate  (micro-stra in/min ) 
vs  time  (min)  Test  UT6F2 
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Axial  Stress  =  2.39  MPa 


Figure  B.9  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  UT6D 
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Axial  Stress  =  2.39  MPa 


Figure  B.10  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  UT6D2 
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LOG  TIME 


Axial  Stress  =  2.92  MPa 


Figure  B.11  Logarithm  plot  of  strain  rate  (micro-st ra in/min ) 
vs  time  (min)  Test  UT6D3 
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Axial  Stress  =  1.86  MPa 


Figure  B.12  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A 
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Axial  Stress  =  2.66  MPa 


Figure  B.13  Logarithm  plot  of  strain  rate  (micro-strain/min ) 
vs  time  (min)  Test  CT6A3 
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